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Preface 


This  report  consists  of  three  papers  describing  research  on  unsteady 
transonic  flow  theory.  The  first,  entitled  "Unsteady  Local  Linearization 
Solution  for  Pulsating  Bodies  at  = 1",  by  S.  S.  Stahara  and  J.  R. 

Spreiter  has  been  published  in  the  AIAA  Journal,  vol.  14,  July  1976, 
pp.  990-992.  The  second,  entitled  "Unsteady  Local  Linearization  Solution 
for  Pitching  Bodies  of  Revolution  at  M = 1:  Stability  Derivative  Analy- 

sis", by  S.  S.  Stahara  and  J.  R.  Spreiter  has  been  published  in  the  AIAA 
Journal,  vol.  14,  October  1976,  pp.  1402-1408.  The  third  publication  is 
ICAS  Paper  No.  76-06,  entitled  "Developments  in  Transonic  Steady  and 
Unsteady  Flow  Theory",  by  J.  R.  Spreiter  and  S.  S.  Stahara.  This  paper 
was  presented  by  Professor  J.  R.  Spreiter  at  the  10th  Congress  of  the 
International  Council  of  the  Aeronautical  Sciences  (ICAS)  in  Ottawa, 
Canada,  October  3-8,  1976. 

These  papers  are  based  on  research  carried  out  under  Contract  No. 
N00014-73-C-0379,  Project  No.  NR  061-215  sponsored  by  the  Office  of  Naval 
Research  with  Mr.  Morton  Cooper  as  Technical  Monitor. 
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Unsteady  Local  Linearization  Solution 
for  Pulsating  Bodies  at  Mx  = 1 

Stephen  S.  Stahara' 

Nielsen  Engineering  <£  Research,  Inc. , 

Mountain  View,  Calif. 
and 

John  R.  Spreitert 

Stanford  University,  Stanford,  Calif. 

Introduction 

WITH  the  continuing  development  of  successful 
techniques  for  solving  steady  transonic  flows,  con- 
siderable interest  has  been  focused  recently  on  the  develop- 
ment of  methods  to  solve  unsteady  transonic  problems.1  In 
this  note  we  describe  the  local  linearization  solution  for  tran- 
sonic flow  past  slender  bodies  of  revolution  undergoing 
oscillatory  pulsatile  motion  of  the  body  surface.  This  result 
provides  the  fundamental  unsteady  source  solution  from 
which  higher-order  multipole  solutions  (dipole,  etc.) 
necessary  to  describe  more  complex  unsteady  motions  (e.g., 
translation,  rotation)  can  be  obtained.  The  theory  is  based  on 
the  concept  of  dividing  the  flow  into  steady  and  unsteady 
components  and  solving  the  resultant  equations  by  the  local 
linearization  method. : The  analysis  is  developed  generally  for 
sonic  and  near  sonic  flows,  with  specific  applications  made  to 
parabolic  arc  half-bodies  and  cones  at  freestream  Mach  num- 
ber M „ = 1 . The  results  indicate  the  correct  convergence  to 
nonlinear  quasisteady  theory  as  the  reduced  frequency  of 
oscillation  based  on  body  length,  k— 0,  and  to  linear  acoustic 
theory  as  k becomes  large  (k  2 I ).  For  k <1,  a range  of  prime 
importance  in  many  flutter  and  stability  applications,  the  un- 
steady solutions  exhibit  a significant  nonlinear  thickness  ef- 
fect induced  by  the  steady-state  solution,  much  like  that 
displayed  in  the  two-dimensional  case.  ’ This  indicates  a basic 
shortcoming  of  linear  theory  in  this  frequency  range. 

Analysis 

The  concept  that  a major  body  of  transonic  flow  problems 
can  be  predicted  accurately  within  the  framework  of  inviscid, 
nonlinear  small-disturbance  theory,  described  by  the  equation 

, I I 

( / — Af  ® ) 0, , + 1 1>„  -t-  — 0,  + , 0ttn 

r r 

= Mi(y  + 1)0,0,,  +\fiD0n  + 2M„0,  <1) 

has  been  well-established.4'  In  Equation  (I),  is  the 
freestream  Mach  number;  (x.r.S)  are  nondimensional  body- 
fixed  cylindrical  coordinates  with  (v.M  normalized  by  body 
length  I’and  with  the  x-axis  directed  rearward  and  aligned  with 
the  body  centerline;  / is  nondimensional  time  normalized  by 
t/U* , where  V F is  the  freestream  velocity;  > is  the  ratio  of 
specific  heats  equal  to  7/5  for  air;  and  e>  is  the  dimensionless 
perturbation  velocity  potential  Although  a variety  of  sub- 
cases of  Eq.  (I)  exist4"  depending  upon  whether  the  time 
behavior  of  the  motion  is  very  slow  (quasisteady:  0,  , 0, 
neglected),  somewhat  more  rapid  (mildly  utnsteady:  0 
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neglected),  or  very  rapid  (linear  high  frequency:  0,0, 
neglected,  the  authors  have  considered  the  more  general  case, 
encompassing  all  frequencies  given  by  Eq.  (I ). 

For  the  general  axisymmetric  oscillatory  flows  conquered 
here,  it  is  convenient  to  expand  the  solution  into  a steady  and 
unsteady  component  Thus,  set 

0ix,r,t,t)  =01  (jx.r)  + R.P.(i(jr.r)e‘‘  | (2) 

where  0,.  is  ihe  axisymmetric  steady  perturbation  potential, 
which  satisfies  Eq.  ( I ) with  the  W and  r terms  omitted,  0 is  the 
complex  amplitude  of  the  oscillatory  perturbation  velocity 
potential,  k is  the  reduced  frequency  defined  by  k = u>t'U„, 
and  R.P.  signifies  the  real  part  of  a complex  quantity.  Based 
on  the  assumption  that  small-amplitude  oscillations  are  ap- 
propriate for  flutter  and  stability  analysis,  the  equation  for  0 
becomes 

0,,  + (//r)0,  = (M;-/  + M-;(y  + />*„]*,, 

+ \Miiy  + l)0,„+2iMik)0, -M:,pu0  (3) 

which,  although  linear,  nevertheless,  remains  formidable 
because  of  the  variable  coefficients  and  mixed  elliptic- 
hyperbolic  type 

The  boundary  condition  at  the  body  surface  can  be  decom- 
posed analogously.  Upon  setting 

Rix.t)  =eR(.v)  +R.P.|6R(x)e,4'J  (4) 

where  (R,  R)  are  normalized  functions  describing  the  steady 
and  oscillatory  components  of  the  body  ordinates,  and  (e,5) 
are,  respectively,  the  normalized  maximum  body  thickness 
and  the  dimensionless  amplitude  of  the  unsteady  oscillations, 
one  obtains 

0h  ix.R  1)  = Ri  + Oie  't  nt)  (5a) 

0,  (jt.R,)  =i£  R‘  + ikR  + -S—  R J *-0(5t  t’ nr6  ) (5b) 

where  R/  = «R.  and  primes  indicate  differentiation  with 
respect  to  x.  The  first  two  terms  on  the  right-hand  side  of  Eq. 
(5b)  are  familiar,  since  they  also  appear  in  the  oscillatory  thin- 
wing  problem.5  For  the  slender  body  case,  however,  the  third 
term  arises  from  the  imposition  of  the  no-flow  boundary  con- 
dition at  the  actual  oscillating  body  surface. 1 where  a Taylor 
series  expansion  about  the  mean  position  r=R,  is  used  to 
remove  the  resulting  implicit  dependence  on  5.  Finally,  the 
corresponding  expressions  for  the  surface  pressure  coef- 
ficients are 

CPl{x,R/)=-20„{x,Rl)-R,1  (6a) 

Cp(x,R,)  = -.?[<»,  (Jr ,R,)+ik0{x.R,) 

+ (r;+^-)  r + r;<r-  +/*R)]  (6b) 

s R 1 t J 

Use  of  the  method  of  matched  asymptotic  expansions  serves 
to  identify  the  logarithmic  behavior  of  both  the  steady  and 
unsteady  components  near  the  body  axis.  Consequently,  the 
differential  equation  and  surface  boundary  condition  for  the 
unsteady  component  can  be  expressed  in  the  compact  form 

\,0„  +\;0,  +\,0=0„+  (Urt0,  (7) 

lim  (x0,  )=g(Jr)=f6R[R'+ikR+<R'/R)R]  (8) 

r -0 

where  (\,.  \2.  k, ) can  be  identified  from  Eq.  (3). 
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Three  fundamentally  different  differential  equations  and 
solutions  occur  for  <£,  depending  upon  the  sign  of  A(;  i.e., 
whether  \,<0  (subsonic),  A,  = 0,  (sonic),  or  \,>0  (super 
sonic).  The  regions  are  illustrated  in  the  upper  left  of  Fig  I 
for  the  forepart  of  a slender  convex  body  at  M , = 1 . The 
solution  for  the  sonic  region  2 must  merge  continuously  with 
that  for  the  subsonic  region  I ahead  of  it  and  the  supersonic 
region  3 behind  it. 

The  procedure  for  determining  the  local  linearization 
solution  for  the  unsteady  component  © requires  first  that  the 
steady-state  solution  <f>,  be  obtained  to  evaluate  the  variable 
coefficients  (A,.  A,).  Then  an  approximate  solution  for  © is 
determined  by  replacing  the  variable  coefficients  (A,,  A..)  tem- 
porarily in  Eq.  (7)  by  constants,  and  solving  the  simplified 
equations  that  result  in  each  of  the  three  regions  identified 
previously.  In  the  original  procedure  for  steady  transonic 
flow,  the  next  step  would  be  to  calculate  the  surface  ac- 
celeration d:4>[x.R  i )/dx: , replace  the  constants  (A,,  A.)  by 


oix.r) 

y.ipci  sonic 


: ( [*(*)•  j£„,(.-M)  +£„(/»  Jr)  +f  n(^r  )] 

*1 <^gU)[e-^  +e  «**-*•)«] 

(A  <>0)  (II) 


where  £,„(Z)  is  the  complete  exponential  function  of  complex 
argument  Z.  * C is  Euler's  constant, 

rl„  = Ai/A.„  A,=  [Aj  ]/2A, 

and 

A;  = [A.  + (A;!-<A,Aj  ) ]/2A, 


fig.  1 I'nsleads  pressure  distributions  on  various  parab»ltc-an  hatf- 
hodies  undergoing  pulsatile  surface  oscillations;  local  lineari/alion 

, quasi  steads  — , acoustic  theory  A -O.OI.  • . f-d.f. 

• ; A = 1.0,  a;A=2.(l.  • ;A=5.I)« 


the  functions  that  they  original  represented,  and  finally  in- 
tegrate the  resultant  second-order  ordinary  differential 
equation  along  the  body  surface  to  obtain  both  diix.R  ) dx 
and  ©(.x, R ,).  to  be  used  in  Eq.  (6b).  In  the  results  reported 
here,  the  authors  have  used  a variant  of  that  procedure  by 
calculating  ©,  (x.R ,)  and  <^(x.A?,)  directly  from  the  simplified 
equations.  Near  the  body  surface,  these  solutions  are 


r 

/A,r2\  ) 

[*<*)• 

£in(zl,sr)+/n(~)+rj 

t"g(*)  ~ g(( ) 
+ )o  x-( 

4>(x.r)  = [ [g(x)« 

Mlbsilllii  2 l 

. ( ~X,r' 

+ fnl  

\ 4x0 -x) 

t',g({)  »|W 

■J.  {-v 


,-tid{J  (X,  = 0) 
[e,„M,x)  +E„l-A,U-x)] 
« 

e4z«t  “df]  (A ,<0) 


<9) 


(101 


Results  determined  by  employing  the  solutions  given  by 
Eqs.  (9-11)  in  their  respective  domains,  are  shown  in  Fig.  I. 
Exhibited  in  the  two  lower  plots  are  the  normalized  magnitude 
and  phase  (in  degrees)  of  the  unsteady  surface  pressure 
distributions  Cp  for  a parabolic-arc  half-body,  with  <=0.10 
executing  pulsatile  oscillations  of  its  body  surface  propor- 
tional to  the  local  radius  (R  = R).  In  those  results,  as  well  as 
others  presented  here,  the  steady  solutions  required  as  input  to 
the  unsteady  calculation  were  determined  by  the  local 
linearization  method  for  axisymmetric  bodies,2  which  is 
known  to  provide  good  accuracy  for  the  shapes  being  con- 
sidered. Also  indicated  on  those  plots  are  the  results  provided 
by  quasisteady  theory  and  by  linear  acoustic  theory,  to  which 
the  present  results  converge  for  small  and  large  k.  The  close 
correspondence  between  the  nonlinear  and  quasisteady  results 
for  k =0.1  implies  a substantial  nonlinear  thickness  effect  of 
the  steady  flow  upon  the  unsteady  component  at  low  frequen- 
cies, and  indicates,  furthermore,  the  quasisteady  theory  can 
provide  good  results  in  this  range.  In  contrast,  the  com- 
parisons with  acoustic  theory  indicate  large  discrepancies  for 
small  k.  particularly  in  phase  angle,  which  tend  to  disappear 
only  when  k is  approximately  2. 

A further  evaluation  of  the  low-frequency  results  predicted 
by  the  present  method  is  provided  by  the  plot  in  the  upper 
right  of  Fig  I The  pressure  distributions  indicated  by  the  un 
steady  analysis  for  a basic  parabolic -arc  half-body,  having  a 
maximum  thickness  ( 0.10,  undergoing  slow  expansion  to  r 

0.12  and  contraction  to  < = 0.08,  are  compared  with  results 
predicted  by  the  local  lineati/ation  theory  for  steady  flow  past 
three  sush  bodies.  The  results  are  in  good  agreement,  in  spue 
of  the  linearization  of  the  unsteady  component  and  the  sub- 
stantially different  methods  of  solution 

Figure  2 shows  the  analogous  results  for  a cone  with  i 
= 0.10  Although  the  analysis  derived  here  applies  strictly  to 
smoothly  accelerating  flows  on  continuous  bsidies  a con- 
dition severely  strained  for  the  flow  in  the  vicinity  of  the 
shoulder  on  a cone  application  ol  the  previous  techniques  is 
still  possible  by  adopting  the  strategy  of  using  the  unsteady 
sonic  solution  predicted  by  the  current  method  along  the  cn 
tire  length  of  cone  This  is  plausible,  since  the  subsonic 
solution,  which  normally  would  be  joined  to  the  sonic  result 
for  points  in  the  region  ahead  of  the  shoulder,  in  last  becomes 
ill-conditioned,  this  is  because,  for  a cone,  the  sonic  point  is 
fixed  at  the  end  of  the  bods  Again,  the  quasisteady  predic 
lions  of  the  present  method  shown  in  the  upper  right  hand 
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tig.  2 l nsteady  pressure  distributions  on  various  cones  undergoing 
pulsatile  surface  oscillations;  local  linearisation  . quasi-stead) 

, acoustic  theory  k = 0.1,  • ; A = 0.5.  ■ ; A = 1 .0,  * ; A = 2.0, 

A - 5.0.  ♦. 

plot  of  Fig.  2 indicate  very  good  agreement  with  the  steady 
state  results,  whereas  the  comparisons  with  the  acoustic 
theory  results  shown  in  the  bottom  plots  indicate  a much  more 
rapid  approach  to  acoustic  theory  as  the  frequency  increases 
for  a cone,  in  contrast  to  a parabolic-arc  half  body. 
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Unsteady  Local  Linearization  Solution  for  Pitching 
Bodies  of  Revolution  at  M =1:  Stability  Derivative  Analysis 


Stephen  S.  Stahara* 

Xielsen  Engineering  & Research,  Inc.,  Mountain  View,  Calif. 


John  R.  Spreitert 

Stanford  University,  Stanford,  Calif. 


An  account  is  provided  of  recent  theoretical  results  lor  unsteadv  inviscid  transonic  flows  about  axisvmmelric 
bodies.  I he  unsteady  local  linearization  solution  is  developed  to  calculate  the  unsteady  pressure  distributions, 
for  sonic  or  near-sonic  freest  ream  flows,  on  the  surface  ot  several  slender,  pointed,  axisvmmelric  bodies  un- 
dei„.ing  low -frequency  oscillatory  pitching  motion.  Kesults  determined  for  static  and  dynamic  stability 
derivatives  exhibit  good  agreement  with  the  limited  experimental  data  available. 


Introduction 

THE  purpose  of  this  paper  is  to  describe  a theoretical 
procedure  for  determining  the  unsteady  flovvfield  at 
freestream  Mach  numbers  at  or  near  one  about  slender 
axisymmetric  bodies  undergoing  oscillatory  pitching  motion. 
The  analysis  is  based  on  the  small-disturbance  theory  of  in- 
viscid unsteady  transonic  flow  and  makes  use  of  the  concept 
of  splitting  the  flow  into  steady  and  unsteady  components  and 
solving  the  resultant  unsteady  potential  equation  by  the  local 
linearization  technique.  An  expansion  of  the  unsteady 
solution  is  obtained  for  low  reduced  frequencies,  for  which 
the  nonlinear  thickness  and  Mach  number  effects  of  the 
steady  flow  are  significant,  and  the  Adam-Sears  iteration 
procedure  is  used  to  determine  the  unsteady  near-field  poten- 
tial. from  which  flow  properties  on  the  body  surface  are 
calculated.  The  approach  generalizes  the  results  of  Landahl, 1 
Liu,:  Liu  et  al.. ' and  Ruo  and  Liu, 4 both  in  the  sense  of  en- 
compassing the  entire  low-frequency  range,  including  A =0. 
and  also  in  eliminating  the  “parabolic”  assumption/4 
thereby  properly  accounting  for  local  regions  of  subsonic, 
sonic,  and  supersonic  flow.  Applications  are  made  to  various 
parabolic-arc  half-bodies  and  cones  at  A/,  = 1 , and  results  are 
presented  for  unsteady  force  and  moment  distributions  and 
static  and  dynamic  stability  derivatives.  The  examples  were 
selected,  insofar  as  possible,  to  enable  comparison  with 
existing  data  and  other  theories.  We  note  that  alternative  ap- 
proaches, employing  the  recently  developed  methods  of 
lsogai‘  and  Dowell/  that  arc  based  on  the  same  general 
notions  of  the  local  linearization  method  are  also  possible  and 
would  provide  an  interesting  further  comparison  of  these  dif- 
ferent procedures. 

Analysis 

Basic  I-  qualions 

The  analysis  is  developed  in  terms  of  a body  fixed  cylin- 
drical coordinate  system  centered  at  the  nose  with  the  x axis 
directed  rearward  and  aligned  w ith  the  longitudinal  axis  of  the 
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body,  and  the  0 = 0 deg  axis  directed  to  the  right,  facing  for- 
ward. A slender,  pointed  body  of  revolution  immersed  in  a 
zero  angle-of-attack  steady  flow  is  assumed  to  undergo  small- 
amplitude  harmonic  pitching  oscillations  about  a point 
located  on  the  body  axis  at  x = a.  With  the  fundamental 
assumption  that  the  body  is  sufficiently  smooth  and  slender  so 
that  the  inviscid  small  disturbance  theory  applies,  a per- 
turbation velocity  potential  o may  be  defined  according  to' 

o (x.r.O.t)  = U*[  (x-a)  cos 6 + r sin0  sin6  + <t>(x,r,0,t)  1 


where  Ux  represents  the  freestream  velocity.  The  angular 
displacement  6 (/)  is  given  by 

5(f)  - RP[bf,e‘k'  ] (2) 

where  6(,  is  the  maximum  displacement,  RP  signifies  the  real 
part  of  a complex  quantity , k is  the  reduced  frequency  defined 
by  k = u)l’/Ux , and  o satisfies 

I I 

(/  — ) <2> , , +£,,  + 0,  + . = M i ( > -f  / )0, O, , 

r r 

+ , +2M:x<t>xl  (3) 

In  Eq.  (3).  \lx  is  the  freestream  Mach  number;  (x.r,$)  are 
nondimensional  body-fixed  cylindrical  coordinates  with(.r.r) 
normalized  by  body  length  ft  t is  nondimensional  time  normal 
ized  by  i/Ux  ; and  > is  the  ratio  of  specific  heats  equal  to  7 5 
for  air. 

For  the  oscillatory  flows  considered  here,  it  is  convenient  to 
expand  the  solution  into  a steady  and  unsteady  component. 
Thus,  we  set 

6(x.r,6.t)  = (*,r)  + RP\4>{x.r)  sinfle'1']  (4) 

where  o is  the  axisymmetric  steady  perturbation  potential, 
which  satisfies  Eq  (3)  with  the  (0.0  terms  omitted,  and  o is 
the  complex  amplitude  of  the  oscillatory  perturbation  velocity 
potential.  On  the  assumption  of  small  amplitude  oscillations 
appropriate  for  flutter  and  stabilitv  analysts,  the  equation  for 
o becomes 

<j>  + 1 O. — . 0=  I M v,  - / + Mi  (>  4-  /)$  . |o. 

r r * 

♦ V#-:  <>  + /)*,  -Uik-6  (5) 


mscoi®  rmEB 


which,  allhough  being  linear , is  formidable  lo  solve  because  of  the  variable  coefficients  and  mixed  elliptic -hyperbolic  t)  pc. 

The  boundary  condition  of  flow  tangency  at  the  body  surface  can  be  decomposed  analogously.  In  the  body-rixed  coordinate 
system,  these  conditions  are’ 


<h  (X.R)  = R +0(t  > Inr)  (6a) 

<*>,  (x.R)  =S„  [/+*(* -a)  I +R  i,  U.R)  -/W„RR  -t- 0(&„* J Inf.ij)  (6b) 

where  R = rR(x)  describes  the  body  ordinates.  Here,  «■  is  the  normalized  maximum  body  thickness,  and  R is  a function  of  order 
one.  At  infinity,  the  perturbation  potentials  are  required  lo  vanish  in  an  appropriate  manner  Since  attention  will  be  confined  to 
flows  with  M = I having  all  shock  waves  downstream  of  the  region  for  which  results  will  be  calculated,  we  will  not  require  the 
corresponding  relation  for  conditions  on  the  two  sides  of  a shock  wave.  The  preceding  equations  provide,  therefore,  the  fun- 
damental relations  for  the  analysis  to  follow. 

Local  l inearization  Solution 

The  differential  equation  (5)  for  the  unsteady  component  can  be  expressed  in  the  compact  form 

- 

A,< j>, , + Aj0  = 0„  + — <t>, j 


\,  = Mi  -I  +Mi(y  + /)$/,.  ^2  = MiiH  + +i2kMi,  X3  = -k2Ml 

Three  fundamentally  different  differential  equations  and  solutions  occur  for  © depending  upon  whether  the  coefficient  X,.  which 
corresponds  in  the  small -disturbance  approximation  to  - I.  is  negative,  near  zero,  or  positive.  * The  zone  of  influence  of 

each  solution  is  different  for  each  case,  as  is  appropriate  since  they  correspond  to  local  regions  of  subsonic,  sonic,  or  supersonic 
flow.  Each  solution  then  forms  the  basis  of  the  ensuing  analysis  in  its  appropriate  region,  with  the  requirement  that  each  partial 
solution  is  required  to  merge  continuously  with  that  for  the  adjacent  region.  The  local  linearization  method  proceeds  by  replacing 
temporarily  the  variable  coefficients  (X/tX, ) in  Eq.  (7)  by  constants  and  solving  the  simplified  equations  that  result  in  each  of  the 
three  regions  discussed  above.  Next,  the  quantities  ©.  (x,R)  and  4(x.R),  which  are  required  for  evaluating  surface  properties  of 
the  How,  such  as  the  unsteady  pressure  coefficient,  are  determined  by  first  obtaining  them  from  the  simplified  equations,  then 
replacing  the  constants  (X/(X*)  by  the  functions  they  originally  represented,  and  finally  evaluating  the  resulting  expressions  at 
points  along  the  body. 

The  three  starting  solutions  for  0 are  given  by 


g(S)  exp  - 


di  (X,=0) 


/ d C1 
<t>(x,r)  = - — l 

MltMHIk  * or 


/ a f'  s(i)  expl  -A,  (x-j)  +/4_'|  (x— j)  •’  + ( —6,  )f’  I’  ’ 1 


|x-J)  •’  + ( -X,  )r’] 


di  (X,<0) 


*(;r,r)  = 7 1 f R(Oexp1-^^-f)Jco.h|^Md(Jr-Q^-A^  ]^  (X|>0)  (10) 

2 dr  l(x  — {'  — X,r- 1 ' 

where  g(x)  is  the  doublet  strength  to  be  determined  later  from  the  boundary  condition, 

/l0  = Xi/X_,,  /l,  = Xj/2X,,  and  A;  =x/X| -4X,X,  /2\, 

In  determining  the  local  linearization  solution,  the  forms  of  these  solutions  in  the  near  field  are  required.  These  could  be  obtained 
by  expanding  the  above  solutions  directly  for  small  r;  however,  a more  convenient  method  is  to  Fourier  transform  the  basic 
equation  (7)  for  the  three  different  regions  with  respect  to  X.  solve  the  transformed  equations,  expand  the  Kernel  functions  in  the 
transformed  solutions  for  small  r,  and  then  transform  the  simplified  results  back  to  the  physical  plane.  Use  of  either  technique 
leads,  after  much  manipulation,  to  the  following  forms: 


Ole.'-)  = - [***'  + ^ [a (.off,,,  (A,*)  +ln(~4X')  +c-/| 

+ j|*  \(Ile  soi.  t'rf(-X,*(0>?-j-  j]  +0(r'lnr)  (X,=0 1 


■Mx.r)  + ''(A(x)(£, „H,«r)+£, „(-/!, {/-arjl+lnl^^^]-/} 

v.ihvonix  * r 4 

+ j*  AU)-A(I),.  4„.  j'  "d( 

+ ([X'C  a +<*')  -Xz]®<0*  — X/S  + {{X' ( -/Ar)  +X.]*( ')  + X,g'  <0,| }]  +0,r,,nr)  <X,«M  (12) 
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f'  A (x)  -A(£) 

1 e 

x-i 


©(x.r)  = -l  — + ' [a«>(£„  (-4.,  x)  +£„W()  +|n(~r)-']  + L 

'uprr%i>nK  ‘ X *»  ' ' 

„[[*,('  +.4,)  .o  + [x,('  +4,)-X..]e  -+»  ^ ) *^  ]] +0,r ' Inr)  <X,>0>  (13 


where 


AUI  s=X,g"  (x)  + X.,g’  (x)  + X jg (x ) 

£„.  (Z)  is  the  complete  exponential  function  of  complex  argument  Z , C = 0.577216  is  Euler  s constant, 

A,  = |X..  - \X::  -4X;X,)/2X, , 


and 


/4  4 — ( X > + v X % — 4X ; X « ) /2X  • 


Low  Frequenc>  Solution 

Of  primary  importance  in  this  investigation  are  low  frequency  oscillations,  that  is,  k = 0 (<  * ln< ) , where  the  nonlinear  thickness 
and  Mach  number  effects  of  the  steady  flow  strongly  influence  the  unsteady  flow  For  this  frequency  regime,  it  is  convenient  to 
decompose  the  unsteady  component  in  the  follow  ing  manner 

0(x,r)e'k  =60'  (x,r)  + 60  1 (x.r)  (14) 


where  6 = /A:6and  the  superscripts  (i,o)  denote,  respectively,  the  in-phase  and  oui-ot'-phase  unsteady  components. 

What  remains  now  is  to  obtain  to  appropriate  order  the  doublet  distribution  g(x)  appearing  in  Eqs.  (11-13).  An  effective  means 
employed  in  similar  investigations'4  is  to  make  use  of  the  Adam-Sears  technique.  If  we  represent  the  doublet  distribution  and 
potential  by 

g(x)  =g;  (x)  +g„>  (x)  <ls> 

© ' (x.r)  = ©'/"  (x.r)  +0."  (x.r)  U6> 


where  g,  (x)  and  ©'/'  (x.r)  represent  the  first-order  effects,  whereas  g,  (x)  and  ©V'  (x.r)  arc  of  a higher  unspecified  order. ap- 
plication of  the  surface  boundary  condition  equation  (6)  serves  to  determine  the  distributions  (g,.g: ) explicitly,  and,  from  these, 
the  solution  for  the  near-field  potential  follows  directly.  The  result  of  these  operations-after  taking  into  account  that 
X,  =0(('lnr),X.  = 0(<  Tor ) + 0(/t),X,  = 0(k:).6  = ikh,  ignoring  0<*- - > terms,  and  retaining  terms  in  the  near-field  potentials  to 
0(e  '.r  ;ln< ) as  is  consistent  with  the  boundary  condition  equation  (6)  - is  that  g , - 0(t  - ),©  " =0(r - /r)  represent  the  slender -body 
doublet,  whereas  g.  =0b  ‘.i  4 ln<  ),©v  =0(< 4 Ir.i  'r  lnr.<  -r)  represent  the  second-order  correction  due  to  body  thickness,  reduced 
frequency,  and  Mach  number.  Specifically,  we  find  that 


© © ( 


© ( t.r)  ©.(x.r) 


x.r)  Otx.r) 


+ 0(<  4 l’«  ) for  r = Ot<  I 


(17) 


w here 


and 


o ( v.r)  - 


/ 0(0 

4n  I r 


° 1 QiX) 

O X-O  xr 

Q (x) 

©.=  - 

4x  • r 

(t?) 


( T / IX.  I x V X..„,.  ' 1 ■■■  Q {x)-Q  U) 

i )+c  +,+  ,*.v f,x  J+L  x~i 


^ ) 


+ A©  (x.r) 


(18) 


<t9) 


(20) 


0 (x)  -Q •({) 


/ O(x)  / ( f x,  1 , X ,)  I"  Q (x)  - 

© I x.r)  + r 2M  , (q  (x)  (£„  - ( / — x)  + ln  - + 1 j + \ 

4,  nr  ' I X «/'  4xll  — x)  > .'ii  x- 


di 


\ Q l(*  0 e"  ' ‘ “ ■ d{)  +fx.0'(x)+X.  Q (x) 

.’■t  | — x nr  f 


X, 

(/-X) 

X| 


x 10- (£)-<?•  (XH  + X.^iy  (£)  -<?  (X)  | 

e 

X, 

rfc  + 

+o.,„l 

|r  ' ‘ 1 

1 - < 

k 

(/-X) 

[ l-x  V 1 

1 « 1 

+ A©  (x.r) 
(21) 
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&Hx.r)  = ' \ y<A)  +f|  \q'  (xlfo,  (~x)  +ln(^r)+/] 

...pmomc  I * r Jl  I L 'h,  /bp  ' 4x-  ' > 

+ i,  0<^lfU>  |/+<?  lVN,l,('l/{]  + fX'<?'(Jt) 

♦ x,  Q ix)}\  —St*\  I +r  Mg-U)-Q-(»>^iO-  U) -O’ <111^]^^, 

X-  .’1/  Jf-f  J 


- , (x-a)Q  Ux)  MIQ'(X)  I Q(x)  Q(x)  ,1 

AO(x,r)=-- = + - --- 1 In- + r Inc-  (23) 

4x2r  2x  ' xr  x J 

In  Eqs.  (20)  to  (22),  the  coefficients  (X,,X , ) are  evaluated  on  the  body  surface  r = R and,  consequently,  are  functions  of  x only, 

\,KP  = RP(\2)=ML(y  + l)<t>,"  (x.R)  (24) 

Q(.v)  represents  the  body  area  distribution,  and  is  given  by  Q(X)  = xR2  ( x ) and  Z/3 signifies  the  imaginary  part  of  a complex  quan- 
tity. Also,  implicit  in  those  equations  arc  the  appropriate  limiting  forms  as  A — 0.  For  example. 


= I +0(k2 ) 


tan 

lim 

l(2kM- 

) 

* 

/2*A(  . v 

1 ’ 

i 

i >* 

2.W-;  t 

k 

... 

W 1 

We  note  that  the  present  sonic  region  results  are  in  agreement  with  the  low-frequency  parabolic-method  results  ot  Refs.  3 and  4 
when  several  higher  order  terms  (in  O inconsistently  included  in  those  results  are  omitted.  With  4>(x,r)  known,  the  surface 
pressure  coefficient  can  be  determined  for  low-frequency,  small-amplitude  oscillations  from  the  result  that 

C„(x.R.d,t)=CPl  (.*./?)  + | f>C'r(x,R  ) +6C”  (x.R  ) |sin0  (27) 

where  the  steady  pressure  coefficient  is  given  by 

CPl(x,R)=  -2<t>,x  (x.R)  -R  2 (28) 

while  the  in-phase  and  out-of-phase  unsteady  pressures  arc 

0<x.R)  = -2o„,  (v.R)  +2[l  {X.R)-(l-My  )r  -’I*/,  (x.R) -0V,  (x.R ) J = C;,(  (x.R  ) + 0,  (x.R ) (29) 

O (x.R ) = -2|0';(X.R)  VO,  (x.R  ) ) + 2[|(/-M'.  >0  (<  (x.R) -(/ -‘M*  ) R'  •’  |o  ; (x.R ) 

+ IMio,  (X.R)  + 1 MiR'  • |0,  (x.R)  +Ro,  (x.R)  + R R ' -0^  (x.R)  -Oj  (x.R ) ] =C",  (x.R ) +C;:, (x.R ) (30) 

where  (O';  0‘,‘  represent,  respectively,  the  first-  and  second-  order  unsteady  pressure  coefficients.  From  these,  the  flutter  and 
stability  derivatives  can  be  calculated  directly  using 


c„  I cr  f c;,(x,R)  1 f c;,(x,R)  i , 

"*  l = --,  <x-a>  y +R(x)R'(x)  _ ■ dx 

Cw  ''  1 I 0 ?<x.R>  I C/.(v.R>  1 


\ f'A  1 =-  \ R(x)  \ ' ;>  dx 

l oj,..  l c;(x,R) 

here  again  the  subscripts  (1, 2)  serve  to  identify  the  first -and  second-order  components  of  the  various  derivatives. 


C;,(x.R) 
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Fig.  1 Distributions  of  first- and 
i cone  and  a parabolic  half-bod>  v 
i = 0 at  \1X  =1. 


second-order  flutter  derivatives  for 
kith  < =0.1  oscillating  in  pilch  about 


and  parabolic  half -bodies  oscillating  in  pitch  about  a = 0al  V/,  = 1. 


Results 

In  order  (o  demonstrate  the  effects  of  thickness  ratio, 
reduced  frequency,  pilch  axis  location,  and  body  geometry  on 
the  unsteady  motion,  the  results  of  the  previous  analysis  were 
applied  to  various  finite-length  circular  cones  and  parabolic- 
arc  half-bodies  described,  respectively,  by 


R(xt  =ix 
R(x)  =i{2x-x! ) 


0<xsl 


(32) 


Static  and  dynamic  flutter  derivations,  as  well  as  their 
variation  along  the  body  surface,  were  calculated  for  M,  = I. 
In  the  results  presented,  the  steady-state  solulions  required  as 
input  to  the  unsteady  calculations  were  determined  by  the 
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Fig.  3 Pilch  moment  flutter  derivative  for  a 12.5  deg  half-angle  cone 
atVf  =1. 


Fig.  4 Pitch  axis  and  frequency  dependence  of  damping -in -pitch 
normal-force  coefficient  of  a 10  deg  half  -angle  cone  at  = 1. 


local  linearization  method,10  which  is  known  to  provide 
results  of  good  accuracy  for  the  shapes  considered  here. 

Figure  1 exhibits  in  the  left-hand  plots  the  distributions  of 
the  first-  and  second-order  static  (dC^Jdx.dC^jdx) , , and 
dynamic  ( dC / dx.dC „■  /dx) , , flutter  derivatives  for  a cone 
having  a maximum  thickness  e = 0. 1 undergoing  slow 
oscillations  in  pitch  about  its  nose  at  Mx  = 1 . We  note  in  par- 
ticular the  rapid  variation  of  the  second-order  dynamic 
derivatives  (dCsJdx,dC^s/dx) : in  the  vicinity  of  the 
shoulder,  a result  of  the  logarithmic  behavior  of  the  steady - 
state  solution1"  in  that  region.  Analogous  results  for  a 
parabolic-arc  half-body  are  given  in  the  right-hand  plots.  No 
large  gradients  are  present  in  those  results,  but,  as  expected,  a 
smooth  variation  is  indicated  of  the  various  derivatives  along 
the  entire  length  of  the  body.  The  parabolic-arc  results  for  the 
second-order  distributions  of  normal  force  (dC^^/dx) : and 
damping-in-pitch  normal  force  (dC^s/dx),  presented  here 
agree  in  trend  but  differ  somewhat  in  magnitude  with  those 
determined  by  the  “parabolic”  method. 5 Those  results  would 
correspond  to  using  the  sonic  results  of  the  present  method 
over  the  entire  bods  with  the  additional  provision  that  the  sur- 
face acceleration  (see  F.q.  (24)]  is  held  constant  at  some 
representative  value. ,M 

Figure  2 provides  the  thickness-ratio  dependence  of  the 
present  theory  for  the  static  (C\A,C\fs)  and  dynamic  (C\A, 
Cm?,)  flutter  derivatives  for  both  cones  and  parabolic-arc  half- 
bodies oscillating  in  pitch  about  the  nose  at  the  reduced 
frequency  Hr  = 0.03  and  Vf,  = I.  Also  shown  for  comparison 
arc  various  available  theoretical  and  experimental  results. 
These  include  the  results  of  Landahl.1  the  linear  transonic 
method  of  Liu,:  the  “parabolic”  method  of  Ruo  and  Liu,4 
and  the  results  of  slender  body  theory,  together  with  ex- 
perimental data  ot  Wehrcnd"  for  the  normal-force  coef- 
ficient slope  CNa  of  a 12.5  deg  half-angle  conc-with  a flat  base. 
Considering  first  the  static  derivative  comparisons  shown  in 
the  two  upper  plots,  we  note  that  l.andahl's  results  for 
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iCNt.CMi)  are  only  first  order  in  thickness  ratio  and  con- 
sequently identical  to  slender -body  theory.  The  linear  tran- 
sonic theory  results  of  Liu 2 for  (C\a,CWa)  are  identical  to  the 
present  results  for  zero  frequency;  however,  those  results  con- 
tain in  addition  terms  linear  in  k that,  for  the  frequency 
regime  considered  here,  are  of  higher  order  and  do  not  enter 
into  the  present  results.  Comparison  with  results  of  the 
“parabolic”  method  of  Ruo  and  Liu4  for  C,VA  indicate  con- 
siderable discrepancy  from  the  present  results  as  the  thickness 
ratio  increases  and  presumably  is  due  to  the  inconsistent  in- 
clusion of  higher  order  terms  0(< 6 Int/r.rt4  ln«)  in  that 
calculation.  For  the  dynamic  derivative  comparisons  shown  in 
the  two  lower  plots,  large  discrepancies  and  different  trends 
exist  among  the  various  theories.  Because  l.andahl’s  results 
are  based  upon  a lower  order  surface  boundary  condition,  as 
well  as  a less  accurate  force  and  moment  determination,  those 
results  can  be  expected  to  deviate  from  the  more  exact  treat- 
ment used  in  the  present  analysis  and  those  of  Liu2  and  Ruo 
and  Liu4  as  the  thickness  ratio  increases.  With  regard  to  the 
comparisons  of  the  linear  transonic  theory  results  of  Liu,2  we 
note  that,  for  the  parabolic-arc  halt-bodies,  the  failure  of  that 
theory  to  include  the  effect  of  the  nonlinear  steady  flow  on  the 
unsteady  problem  is  sufficient  to  change  the  entire  trend  of 
the  dependence  of  the  dynamic  stability  derivatives  on 
thickness  ratio.  This  provides  yet  another  example  of  both  the 
sensitivity  and  dependence  of  the  unsteady  motion  on  the 
nonlinear  steady  flow  at  transonic  speeds  and  points  out  the 
necessity  of  including  this  effect. 

A final  comparison  of  static  derivative  results  of  the  present 
method  with  experiment  is  provided  in  Fig.  3.  Flere,  the  data 
obtained  by  Wehrend"  for  the  slope  of  the  pitching  moment 
coefficient  CX,A  is  shown  for  a 12.5  deg  half-angle  cone  for 
three  different  pitch  axes  locations  at,  respectively,  63  per- 
cent, 70  percent,  and  77  percent  of  the  body  length  from  the 
nose  and  for  Freestream  Mach  numbers  from  0.65  to  1.30. 
Also  provided  are  the  theoretical  results  of  the  present  method 
for  = I . as  well  as  those  of  Landahl,1  Ruo  and  Liu,4  and 
slender-body  theory.  We  note  again  that  the  linear  transonic 
theory  results  of  Liu2  for  the  static  derivatives  are  essentially 
identical  to  the  present  theory  and  have  been  omitted  for 
clarity  of  presentation.  Good  agreement  of  the  present  results 
is  noted  with  both  experiment  and  the  “parabolic”  results  of 
Ruo  and  Liu.4 

Figure  4 displays  the  pitch  axis  and  frequency  dependence 
predicted  by  the  present  method  for  the  damping-in-pitch  nor- 
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mal-force  coefficient  CV/>  of  a 10  deg  half-angle  cone 
oscillating  in  pitch  at  Mx  = 1.  Shown  for  comparison  are  the 
results  of  Landahl, 1 the  linear  transonic  method  of  Liu,*1  the 
“parabolic”  method  used  by  Ruo  and  Liu,4  and  slender-body 
theory.  Significant  differences  among  the  various  results  are 
observed.  Moreover,  the  well-known  logarithmic  singularity 
predicted  both  by  Landahl’s  method1  and  linear  transonic 
theory2  as  the  reduced  frequency  Ar  — 0 is  quite  evident,  and  is 
in  distinct  contrast  to  the  smooth  variation  with  k of  the 
present  theory.  Moreover,  we  note  that  the  results  of  Ruo  and 
Liu,4  which  are  valid  only  for  very  low  (quaststeadyl  frequen- 
cies and  which  do  not  contain  an  explicit  dependence  on  k, 
compare  favorably  with  the  present  theory  near  k = 0. 

A further  comparison  of  the  results  for  the  damping-in- 
pitch normal-force  coefficient  is  provided  in  Fig.  5,  which 
displays  the  influence  of  body  shape  on  both  the  pitch  axis 
and  frequency  dependence  of  C\A-  Flere,  various  theoretical 
results  for  a cone  are  shown  in  the  two  upper  plots  and  are 
contrasted  in  the  lower  plots  with  the  results  for  a parabolic- 
arc  half-body  having  the  same  normalized  maximum 
thickness  < =0.10.  Greater  disagreement  exists  among  the 
results  of  the  various  theories  for  the  conical  rather  than  the 
parabolic-arc  shape,  a fact  already  observed  in  the  thickness- 
ratio  comparisons  of  (Cv  shown  in  Fig.  2.  This 
is  almost  certainly  a result  of  the  particular  behavior  of  the 
flow  in  the  vicinity  of  the  shoulder  of  conical  bodies,  where 
rapid  variations  in  flow  properties  occur  - an  effect  altogether 
absent  in  the  smoothly  continuous  flow  about  parabolic-arc 
half-bodies. 

Results  analogous  to  those  given  in  Fig.  4 and  5 are  pre- 
sented in  Figs.  6 and  7 for  the  damping-in-pitch  moment  coef- 
ficient Cw.  In  Fig.  6,  results  for  the  damping-in-pitch 
moment  coefficient  CVfs  of  a 12.5  deg  half-angle  cone  at 
M,  = 1 arc  exhibited.  In  addition  to  the  theoretical  results  in- 
dicated, experimental  data  of  Wehrend"  for  pitch  axes 
locations  a =0.63,  0.70,  and  0.77  are  also  included  in  the  plot 
on  the  left.  We  note  that  although  the  present  theory  provides 
results  somewhat  closer  to  the  data  than  the  other  methods, 
none  of  the  theoretical  results  can  be  considered  to  be  in  ac- 
ceptable agreement  with  the  data.  A feature  that  introduces 
additional  uncertainty  into  these  comparisons  is  that  the  ac- 
tual flow  past  a cone  with  flat  base  will  contain  a free 
streamline  emanating  from  the  shoulder  whose  predicted 
location  is  beyond  the  scope  of  an  inviscid  theory.  Whereas 
the  direction  of  that  streamline  will  undoubtedly  be  close  to 
freestream  and  its  influence  restricted  essentially  to  the 
vicinity  of  the  cone  shoulder,  inviscid  theory  predicts  that  the 
surface  pressures  will  be  unaffected  by  this  feature  of  the 
flow,  as  the  influence  of  its  effect  is  limited  to  regions  of  the 
flow  field  downstream  of  the  limiting  characteristic 
emanating  from  the  cone  shoulder.  Discrepancies  associated 
with  this  effect,  together  with  other  uncertainties  related  to 
the  dynamic  testing,"  arc  difficult  to  assess  further  without 
additional  experiment  an  4 analysis.  With  regard  to  Fig.  6,  we 
observe  again  the  differences  between  the  various  methods, 
particularly  as  shown  in  the  plot  on  the  right  at  low  frequen- 
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Fig.  7 l)amping-in-pikh  moment  coefficient  for  a cone  and 
parabolic  half-bo d>  with  # =0.1  at  Mx  = I. 


cies  and  note  the  smooth  variation  with  frequency  of  the 
present  results. 

Finally,  an  assessment  of  the  effect  of  body  shape  on  CK1h  is 
provided  in  Fig.  7,  which  compares  the  pilch  axis  and 
frequency  dependence  for  a cone  and  parabolic  half-body 
with  ( =0.10  at  Mx  = 1.  We  note  in  these  results  the  relative 
closeness  of  the  different  methods  for  both  the  cone  and 
parabolic  half -body. 

Concluding  Remarks 

In  this  paper,  we  have  given  an  account  of  a theoretical 
study  to  calculate  unsteady  surface  pressure  distributions  and 
stability  derivatives  for  slender,  pointed,  axisymmetric  bodies 
undergoing  low-frequency  pitching  oscillations  in  a flow  with 
sonic  or  near-sonic  freest  ream  velocity.  The  results  are  shown 
to  be  in  good  agreement  with  the  limited  experimental  data 


available  for  comparison.  Additional  testing  is  needed,  par- 
ticularly for  the  dynamic  derivatives,  to  provide  a definitive 
evaluation  of  the  theoretical  predictions. 

Finally,  we  note  that  the  local  linearization  procedures  by 
which  the  solutions  have  been  obtained  are  not  restricted  to 
the  particular  examples  displayed  in  this  paper  but  possess 
greater  generality.  The  present  results,  in  fact,  provide  yet 
another  application  of  the  method  demonstrating  the 
remarkable  ability  of  (he  technique  to  obtain  both  an  accurate 
and  economical  approximation  to  a difficult  nonlinear 
problem. 
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discussion  of  significant  results,  usually  obtained 
by  more  them  one  method.  We  start  with  two-dimen- 
sional flows , for  which  the  number  and  variety  of 
methods  is  greatest,  and  with  the  Navier-Stokes 
representation  so  that  the  initial  results  presented 
will  provide  a standard  with  which  the  following 
approximate  theories  can  be  judged.  Three- 
dimensional  flows  are  discussed  subsequently. 


Abstract 


Great  strides  have  been  made  in  recent  years 
in  the  analysis  of  transonic  flows  past  wings  and 
bodies , and  the  initial  extensions  to  wing-body 
combinations,  helicopter  rotors,  and  internal  flow 
through  rotating  turbomachinery  have  been  announced 
Although  the  methods  are  too  numerous  and  diverse 
to  permit  detailed  description,  the  salient  fea- 
tures and  results  of  the  more  significant  are 
revived. 


II.  Two-Dimensional  Flow 


Although  the  majority  of  analysis  has  been 
directed  toward  an  inviscid  level  of  approximation 
transonic  flows  are  often  highly  influenced  by 
viscous  phenomena,  particularly  shock-induced 
boundary-layer  separation.  The  most  complete  anal- 
ysis of  such  a flow  is  that  of  Dievert(l2,13)  who 
employs  the  finite  volume  method  to  solve  the  time- 
dependent  , Reynolds-averaged  Navier-Stokes  equa- 
tions which,  written  in  integral  form,  are 


I.  Introduction 


With  the  propsect  of  efficient  transonic  flight 
afforded  by  the  development  of  supercritical  air- 
foils providing  the  motivation,  with  the  results 
'1-5)  of  a previous  era  extending  from  about  15  to 
25  years  ago  providing  the  foundation,  and  with  an 
enhanced  computational  capability  providing  the 
means , enormous  progress  in  transonic  flow  analy- 
sis is  currently  being  made.  During  this  decade, 
several  new  methods  particularly  adapted  for  use 
with  advanced  computers  have  been  developed  to 
calculate  flows  about  increasingly  complex  aero- 
dynamic configurations.  These  techniques  are  not 
limited  to  the  small  disturbance  theory  that  formed 
the  basis  for  most  of  the  earlier  work , but  apply 
also  to  (a)  the  complete  potential  equation  for 
irrotational  flow,  (b)  the  Euler  equations  for 
general  inviscid  flow,  and  (c)  even  the  Navier- 
Stokes  equations  for  viscous  flow.  Some  of  the 
methods,  such  as  the  numerical  time-dependent 
procedures , have  been  applied  to  all  levels  of 
description.  Others,  as  the  hodograph,  are  re- 
stricted to  two-dimensional  steady  flow  governed 
by  either  the  complete  potential  equation  or  its 
small  disturbance  counterpart.  Still  others,  as 
the  integral  equation  method,  have  broader  poten- 
tialities , but  are  presently  developed  primarily 
for  the  small  disturbance  equation.  As  might  be 
expected,  most  methods  have  been  applied  to  the 
small  disturbance  formulation  and  fewest  to  the 
Navier-Stokes  equations.  Under  favorable  circum- 
stances, the  various  methods  yield  similar  results, 
with  the  effects  of  viscosity  being  confined  to 
limited  regions,  so  that  the  results  of  the  small 
disturbance  theory  agree  with  those  of  the  more 
accurate  theories.  Some  of  the  differences,  more- 
over, stem  directly  from  differences  in  computa- 
tional details  such  as  mesh  size,  order  of  accu- 
racy of  the  difference  algorithm,  etc.,  and  could 
be  reduced  at  the  expense  of  greater  computing 
time  or  with  the  development  of  more  effective 
algorithms. 


where 


in  which  p is  density ^ U and  W are  velocity 
components  and  ix  and  iz  are  unit  vectors  parallel 
to  the  x and  z axes,  E = e + (U2  + V2)/2  is  total 
energy  per  unit  volume,  < is  thermal  conductivity, 
VT  is  temperature  gradient,  a and  T are  normal  and 
shear  stresses,  and  n is  a unit  normal  vector. 

In  addition  to  ordinary  viscous  stresses  in  T, 
Diewert  also  includes  turbulent  Reynolds  stresses 
estimated  using  four  different  algebraic  eddy 
viscosity  models.  If  the  viscous  and  Reynolds 
stresses  are  disregard,  Eq.  (l)  reduces  to  the 
integral  form  of  the  Euler  equations  for  inviscid 
flow, 


Because  the  volume  of  material  is  so  large,  and 
reference  can  be  made  to  several  excellent  reviews 
(6-11)  as  wen  as  original  papers,  our  approach 
here  is  to  provide  a broad  overview  by  combining 
brief  accounts  of  the  theoretical  basis  of  the 
various  methods  followed  by  a presentation  and 


The  integral  forms  are  appropriate  for  the 
finite  volume  method  used  by  Diewert  for  viscous 
flow  and  also  by  Rizzifl*4)  for  inviscid  flow;  but 
other  methods  are  usually  derived  from  the 
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corresponding  differential  equations.  The  Euler 
equations,  for  example,  Eire  as  follows  when  written 
in  conservation  form  to  include  the  shock  relations 
directly 


>=  0 (3) 


Since  these  equations  are  difficult  to  solve, 
they  are  often  approximated  by  assuming  the  flow 
to  be  isentropic.  This  is  correct  without  approxi- 
mation for  flows  which  are  continuous,  i.e.,  shock- 
free,  and  a good  approximation  for  flows  with  shock 
waves  if  they  are  not  too  strong.  Under  these 
circumstances , 

(p/P„)/(p/p„)Y  = exp  (( s - soj)/cv]  (U) 

where  y = Cp/cv  is  the  ratio  of  specific  heats 
at  constant  volume  and  pressure,  s is  entropy, 
and  subscript  00  refers  to  free-stream  conditions. 
Then  Eq.  (3)  reduces  to 
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Magnus  and  Yoshihara  } have  introduced  time- 
dependent  numerical  methods  to  transonic  analysis 
by  solving  these  equations. 

A related  assumption,  which  is  also  exact  for 
shock- free  flow  and  a good  approximation  if  the 
shock  waves  are  not  too  strong,  is  that  the  flow 
is  irrotational.  There  then  exists  a velocity 
potential  $ related  to  the  velocity  by 
q = V4>  that  satisfies 
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For  Moo  = Uso/aoc  f 1,  a number  of  alternatives 
exist  for  k,  but  the  widely  used  expression 
k = y + l)l^/Ux>  provides  a good  all-around 
compromise. 

Equations  (l)  through  (7)  must  be  supplemented 
by  boundary  conditions  and  other  relations  to  speci- 
fy both  a mathematical  problem  corresponding  to  a 
physical  application  and  a unique  solution.  Bound- 
ary conditions  must  uescribe  conditicwis  at  infinity 
and  at  the  airfoil  surface.  With  Eq.  (l),  the 
velocity  must  satisfy  the  no-slip  condition  at  the 
airfoil  surface;  whereas  only  the  normal  component 
must  match  that  of  the  airfoil  surface  with  Eqs.(2) 
through  (7).  With  Eq.  (7),  it  is  customary  to 
linearize  the  boundary  condition  and  transfer  it 
to  the  x axis.  Thus, 


(*z>z=0  = “-(ft*  ft) 


for  an  airfoil  having  ordinates  Z(x,t).  Similarly, 
the  exact  Bernoulli  equation  is  used  with  Eq.  (6), 
whereas  an  approximation 
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is  normally  used  with  Eq.  (7).  Equations  (l) 
through  (*>)  do  not  need  supplementary  relations 
for  shock  waves,  but  Eqs.  (6)  and  (7)  do.  To  com- 
plete the  specification  for  a lifting  airfoil,  it 
is  also  necessary  to  impose  the  Kutta  condition  or 
some  counterpart. 

Of  the  many  methods  for  solving  steady  transonic 
flow  problems,  the  hodograph  method  is  unique  in 
that  it  depends  on  a transformation  that  linearizes 
the  governing  equation  without  approximation  by 
interchanging  the  dependent  and  independent  vari- 
ables. In  this  way,  the  steady-state  form  of  Eq. 
(6)  transforms  to 


(a'  -UJ>$i;u-2UW  4,^+  (a2  -W2)^  = 0 (10 

and  its  small  disturbance  counterpart,  Eq.  (7),  to 


where  a = [Uq  - (y  - l)q2/2]*/2  is  the  local 
speed  of  sound,  and  aQ  is  the  stagnation  speed 
of  sound.  When  applied  to  steady  flow  so  that  the 
time  derivatives  vanish,  Eq.  (6)  is  of  elliptic 
type  where  the  local  Mach  number  M = |7|/a  < 1. 
and  hyperbolic  type  where  M = | q. | /a  > 1.  Both  the 
hodograph  and  finite  difference  relaxation  methods 
have  been  used  extensively  to  solve  this  equation. 

Most  transonic  analysis  both  past  and  present 
has,  however,  been  based  on  a small  disturbance 
approximation  to  Eq.  (6).  For  flow  with  free- 
stream  velocity  Too  = £xUa>  Past  a thin  airfoil 
aligned  approximately  with  the  x axis,  a pertur- 
bation velocity  potential  <j>  = 4>  - u„x  esm  be 
introduced  and  shown  to  satisfy 
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in  which  4>  is  the  Legendre  potential  related  to 
the  coordinates  by  pu  = x and  <K,  = z.  The  Jacobian 
of  the  transformation , J = u£(m£  - 1 + ku)  - u|  in 
the  case  of  Eq.  (ll),  cannot  vanish  where  the  flow 
is  subsonic,  but  may  vanish  where  it  is  supersonic, 
thereby  signifying  an  inconvenient  multiple  mapping 
in  the  hodograph  plane  of  a single  point  in  the 
physical  plane.  As  a result,  use  of  this  method 
is  usually  confined  to  the  subsonic  region  and  a 
limited  part  of  the  supersonic  region  unless  the 
latter  is  small  as  at  slightly  supercritical  Mach 
numbers.  The  remainder  of  the  solution  is  calcu- 
lated by  another  method,  such  as  characteristics 
or  finite  differences.  The  hodograph  method  was 
used  extensively  in  the  early  development  of  tran- 
sonic flow  theory (■? ' , but  difficulties  in  imposing 
boundary  conditions  for  an  airfoil  of  specified 
shape  restricted  its  usefulness.  It  has  come  into 
a new  era  of  significance  in  the  design  of  shock- 
free  airfoils  for  which  it  is  particularly  well 
suited. 
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Two  methods  of  solving  the  hodograph  equations 
are  currently  predominant.  That  of  Nieuwland^^) , 
Takahashi , and  Boerstoel(  18)  relies  on  function 
theory  and  integral  transforms;  that  of  Bauer, 
Garabedian,  and  Korn (19 *20 ) is  based  on  analytical 
continuation  into  the  complex  domain  so  that  the 
solution  can  be  obtained  by  solving  initial  value 
problems.  Both  methods  culminate  with  extensive 
numerical  computations. 

The  first  method  to  be  developed  for  shocked 
transonic  flow  that  depends  intrinsically  on  the 
modern  electronic  computer  is  the  time-dependent 
method  of  Magnus  and  Yoshihara(15»21 ,22) . it 
applies  a modified  Lax-Wendroff  difference  scheme 
to  Eq.  (5)  to  compute  from  a knowledge  of  p,  U, 
and  W at  one  time  new  values  for  a slightly  later 
time.  Repetitive  application  provides  a time 
history  of  the  development  of  the  flow  from  an 
assumed  initial  state.  Solutions  for  steady  flow 
are  obtained  by  imposing  steady  boundary  condi- 
tions and  carrying  the  calculations  sufficiently 
forward  in  time  for  the  transients  to  disappear, 
characteristically  a costly  process. 

Relaxation  methods  were  used  in  the  early  tran- 
sonic studies  of  Emmons (23)  in  physical  variables 
and  Vincenti  and  Wagoner (2U ,25 ) in  hodograph  vari- 
ables, but  modern  applications  derive  directly 
from  Murman  and  ColeV^6)  wh0  introduced  the  use  of 
different  difference  expressions  in  the  subsonic 
and  supersonic  regions  in  conformity  with  the 
respective  domains  of  influence.  The  procedure 
was  initially  applied  to  Eq.  (10)  in  which  <f>zz 
is  approximated  by  centered  differences;  and  $x 
and  4)^  are  represented  by  centered  differences 
in  regions  of  subsonic  flow  and  by  one-sided  up- 
wind differences  in  regions  of  supersonic  flow. 

The  resulting  large  set  of  algebraic  equations 
was  solved  numerically  using  a successive  line 
overrelaxation  (SLOR)  procedure. 

It  was  soon  recognized  that  the  use  of  only 
two  types  of  operators  for  <t>x  and  was  insuf- 

ficient since  difficulties  arose  at  the  sonic  line 
and  shock  wave  where  transitions  between  the  two 
forms  had  to  be  made.  The  problem  was  resolved 
by  introducing  a sonic  finite  difference  operator 
corresponding  to  the  centered-difference  expression 
for  4>zz  and  a shock  operator  corresponding  to 
the  sura  of  the  subsonic  and  supersonic  operators. 
The  latter  when  applied  to  the  steady-state  form 
of  Eq.  (7)  rewritten  in  conservation  form 

(1  - M2  - (R/2 )<t>2 ] + ($  ) =0  (12) 

°°  x x y y 


results  in  substantial  improvement  of  the  calcu- 
lated changes  across  a shock  wave (27).  Presently, 
two  types  of  differencing  procedures  are  widely 
used  with  regard  to  the  small  disturbance  equa- 
tion. To  distinguish  them,  those  employing  the 
shock  point  operator  (SPO),  which  assures  conser- 
vation of  mass  at  shock  points,  are  designated 
FCR  (fully  conservative  relaxation);  while  those 
which  do  not  employ  the  SPO  are  termed  NCR  (not 
fully  conservative  relaxation).  Further  compli- 
cations arise  when  the  method  is  applied  to  the 
complete  potential  Eq.  (6)  specialized  to  steady 
flow.  The  centerline  of  the  Mach  forecone  region 
of  dependence  in  supersonic  flow  is  no  longer 
parallel  to  the  x axis,  but  is  in  the  local 
direction  of  the  streamlines.  The  calculations 


fail  to  converge  if  the  differencing  is  done  in 
the  same  way  as  in  the  small  disturbance  theory 
since  the  region  of  dependence  of  the  differencing 
scheme  does  not  always  include  that  of  the  differen- 
tial equation.  However,  Jameson(28)  showed  that 
convergence  could  be  achieved  if  the  finite  diff- 
erence elements  were  rotated  to  allow  for  the 
change  in  direction. 

Improvements  have  also  been  made  in  devising 
more  efficient  algorithms  for  solving  the  large  set 
of  algebraic  equations  provided  by  the  finite  diff- 
erence procedure.  The  SLOR  method  employed  origi- 
nally has  been  supplemented  by  more  efficient 
procedures )^ >8 >9 ,10 ,29 ,30) > and  further  advances 
exploiting  steadily  improving  computer  capabilities 
are  to  be  anticipated. 


The  integral  equation  method,  stemming  from 
Oswatitsch)  31 ) and  Spreiter  and  Alksne(32),  was 
the  first  method  capable  of  providing  results  with 
embedded  shock  waves  for  supercritical  flows . In 
this  method.  Green's  theorem  is  used  to  transform 
the  steady-state  form  of  Eq.  (!)  into  a nonlinear 
integral  equat ion 


u = u.  ♦ vu2 
L 


(x-X)2  - Cz-Z)2 
[(x  - t)2  + (z  - c)2] 2 


d£  d£  (13) 


where  x = x,  z = (l  - Mi)l/,2z,  u = u/ucr  = 
ku/(l-Mi),  and  u^  represents  the  solution  of 
Eq.  (7)  with  the  nonlinear  and  time-derivative 
terms  omitted.  Originally,  the  singularity  at 
the  field  point  where  u is  to  be  evaluated  was 
enclosed  by  an  infinitesimal  rectangle  having  a 
ratio  X of  height  to  width  of  infinity,  in  which 
case  V = 1/2.  Nixon  and  Hancock) 33)  have  subse- 
quently obtained  V = 1/b  by  enclosing  the  same 
point  by  an  infinitesimal  circle.  Ogana  and 
Spreiter' 31* ) have  investigated  this  matter  further 
by  enclosing  the  field  point  in  an  infintesimal 
rectangle  of  arbitrary  X , for  which  v = (l/ir) 
arctan  X.  This  accounts  for  the  earlier  results 
since  V = 1/2  when  Xs00,  and  V = 1/U  when  X = 1. 
These  differences  are  somewhat  illusory,  however, 
since  the  doublet  integral  is  semi -convergent  and 
its  value  depends  on  X in  such  a way  that  the 
sum  of  the  last  two  terms  of  Eq.  (13)  is  indepen- 
dent of  the  shape  of  the  contour. 

Solutions  of  Eq.  (13)  were  sought  initially  by 
introduction  of  a velocity  profile  f of  the  form 
u(x,z)  = u(x,0)f(x,z,'t,=Q,Z)  so  that  the  double 
integral  could  be  approximated  by  a single  integral 
and  solutions  could  be  obtained  with  a modest 
amount  of  hand  computat ion) 32 , 3 1 ' . Although  this 
method  has  substantial  merit)3b),  the  trend(33,37) 
is  to  direct  numerical  solution.  Nixon  and  Hancock 
) 3 3 ) have  shown,  in  addition,  that  a notable 
improvement  for  the  vicinity  of  the  leading  edge 
may  be  achieved  by  using  a modified  u^  obtained 
by  replacing  the  1 undary  condition  of  Eq.  (8) 
with 
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In  addition  to  the  above  methods  that  may  b*» 
refined  system^*  ically  to  pr  vi  ie  an  essentially 
exact  solution  vithin  the  framev  rk  of  small  dis- 
turbance theory,  there  are  several  approximate 
methods  that  can  provide  h urate  res  ii?  , ft  en 
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very  simply.  Principal  among  these  are  the  local 
linearization(38»39'  and  the  parametric  differen- 
tiation'^) methods  Since  they  are  not  in  such  an 
active  stage  of  development,  no  general  description 
will  he  provided;  although  their  results  will  be 
displayed  for  comparison. 

We  turn  now  to  the  results  provided  by  all  of 
the  methods  discussed.  Figure  1 presents  pressure 
distributions  for  an  18-percent  thick  circular-arc 
airfoil  as  measured)^ ),  and  as  calculated  by 
Diewert'  >13)  for  both  inviscid  and  viscous  flow 


Fig.  1.-  Pressure  distributions  for  thick 
circular-arc  airfoil. 

by  application  of  the  finite  volume  method  to 
Eq.  (l).  The  inviscid  solution  agrees  well  with 
experiment  over  the  forward  half  of  the  airfoil, 
but  is  inaccurate  in  predicting  shock  strength 
and  location,  and  the  pressure  level  near  the 
trailing  edge.  When  the  aft  pressure  recovery  is 
strong,  as  for  Reynolds  number.  Re  = 2*L06,  the 
viscous  solution  agrees  well  with  experiment. 

When  the  aft  pressure  recovery  is  weak,  the  theo- 
retical results  disagree,  probably  because  of  in- 
adequate turbulent  modeling  in  the  separated 
region.  At  the  present  time,  these  are  the  most 
complete  comparisons  between  viscous  theory  and 
experiment;  the  differences  should  be  remembered 
in  evaluating  subsequent  comparisons  between 
inviscid  theory  and  experiment. 

Figure  2 presents  a comparison  of  exact  and 
approximate  solutions  for  steady  inviscid  shock- 
free  transonic  flow  past  two  different  supercriti- 
cal airfoils ' 10 ) . It  shows  that  results  obtained 
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furthermore,  that  the  transonic  small  disturbance 
Eq.  (7)<-  ''  can  provide  results  of  outstanding 
accuracy,  a conclusion  long  evident  from  early 
studies v3) , 

(9) 

Figure  3 from  Ballhaus  shows  pressure  distri- 
butions for  tranr onic  flow  with  embedded  shock  wave 
for  an  NACA  6UaU10  airfoil  calculated  by  application 


o — o Small  disturbance 
Euler 


Fig.  3.-  Exact  and  approximate  transonic  pressure 
distributions  for  an  NACA  6bAUlO  airfoil  at 
two  angles  of  attack,  = 0.72. 

(IS  21  22) 

of  the  time-dependent  method'  ’ ’ ; to  the  Euler 

Eqs.  (5)  for  isentropic  flow  and  of  the  finite 
difference  relaxation  method  to  the  small  distur- 
bance Eq.  (7).  Two  versions  of  the  latter  are 
shown  distinguished  by  the  value  of  m in 
k = M[£(y  + 1 )/L)oo.  The  upper  plots  are  for  m = 2, 
in  wide  use  for  many  years(**®K  To  improve  accuracy 
without  complication  Krupp has  proposed  several 
other  values  for  m,  each  of  which  has  certain 
desirable  properties  such  as  the  accuracy  of  the 
approximation  for  the  shock  jumps  or  the  critical 
pressure  coefficient,  Cj.  Among  these  values  is 
m = 1.69  used  in  the  lower  plots  of  Fig.  3.  The 
results  illustrate  the  dilemma  encountered  in  such 
efforts.  The  results  for  m = I.69  are  superior 
for  angle  of  attack,  a = 0,  but  those  for  m = 2 
are  better  for  a = 2°.  To  us , it  seems  preferable 
to  use  m = 2 for  most  applications,  reserving 
other  values  for  restricted  classes  of  applications 
that  emphasize  a particular  feature  of  the  flow 
that  can  be  better  represented  by  another  value 
for  m. 


Fig.  2.-  Theoretical  shock-free  pressure 
distributions  for  two  supercritical  airfoils. 

by  application  of  finite  difference  relaxation 
methods (**2, ** 3)  to  the  complete  potential  Eq.  (6), 
and  the  hodograph  method (** 3, UU)  to  its  counter- 
part Eq.  (10),  are  virtually  identical;  and, 


As  noted  above,  significant  differences  occur 
for  the  vicinity  of  a shock  wave  depending  on 
whether  the  finite  differencing  is  done  with  the 
FCR  or  NCR  procedures.  These  are  illustrated  in 
Fig.  U for  a nonlifting  6-percent  thick  circular- 
arc  airfoil  at  two  different  values  for 
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Fig.  4.-  Supercritical  pressure  distributions  for 
a thin  circular-arc  airfoil  indicated  by  the  FCR 
and  NCR  solutions  of  the  transonic  small  distur- 
bance equations,  by  the  finite  volume  solution 
of  the  Euler  equations,  and  by  experiment. 

Away  from  the  shock  wave,  the  results  are  in  agree- 
ment, but  near  the  shock  wave  they  are  quite  dif- 
ferent. In  particular,  the  NCR  shock  pressure  jump 
does  not  approach  the  theoretical  normal  shock 
jump;  whereas  the  FCR  solution  shows  not  only  that 
the  correct  condition  is  attained  but  that  the 
shock  is  followed  by  a well-defined  reexpansion, 
as  is  proper.  The  FCR  solution  also  indicates  a 
more  downstream  location  for  the  shock  wave  than 
the  NCR  solution.  Superposed  on  the  results  for 
= 0.909  are  experimental  data(48)  and  a finite 
volume  solution(l**)  of  the  Euler  Eqs.  (2).  This 
comparison  confirms  that  the  FCR  shock  location 
agrees  better  with  t.he  exact  inviscid  location 
than  the  NCR  locations;  although  the  latter  agree 
better  with  experiment.  This  paradox  can  be  re- 
solved by  reference  to  Fig.  1 in  which  it  is  illus- 
trated how  viscous  effects  lead  to  a more  upstream 
location  of  an  embedded  shock  than  indicated  by 
inviscid  theory.  In  addition  to  the  differences 
in  the  surface  pressures  illustrated  in  Fig.  4, 
Newman  and  South (49)  have  shown  that  the  stream- 
lines downstream  of  the  shock  are  displaced  out- 
ward substantially  and  erroneously  by  the  NCR 
method  because  of  spurious  mass  addition  at  the 
shock.  The  FCR  method  avoids  this  deficiency  by 
satisfying  mass  conservation  everywhere. 


a time-dependent  method(-50)j  which  is  fully  conser- 
vative. All  three  methods  indicate  essentially  the 
same  pressure  distribution  on  the  airfoil,  but  the 
NCR  method  again  indicates  a location  for  the  shock 
wave,  this  time  detached,  that  is  too  far  forward. 

Although  the  integral  equation  method  has  been 
partially  eclipsed  by  the  holograph  and  relaxation 
methods,  recent  reevaluations (33,36 ,37)  indicate 
that  it  has  considerable  merit.  Figure  6 shows 
pressure  distributions  for  a 6-percent  thick 
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Fig.  6.-  Pressure  distributions  for  a 6-percent 
thick  circular-arc  airfoil  indicated  by  integral 
equation  and  FCR  solutions  of  the  transonic 
small  disturbance  equation. 

circular-arc  airfoil  calculated  using  two  versions 
of  the  integral  equation  method  and  the  small  dis- 
turbance FCR  finite  difference  method.  The  results 
of  Kraft (36)  have  been  determined  using  a velocity 
profile  to  reduce  the  doublet  integral  to  a single 
integral  and  iterating  in  the  manner  of  Spreiter 
and  Alksne(32).  Those  of  Nixon (51)  have  been  cal- 
culated using  his  extended  integral  equation  method 
in  which  the  doublet  integral  is  evaluated  by  divid- 
ing the  region  of  integration  into  a number  of 
streamwise  strips  across  which  interpolation  func- 
tions are  used  to  express  values  for  the  integrand 
in  terms  of  values  along  the  strip  edges.  Similar 
results  for  the  critical  Mach  number  have  also  been 
determined  by  Ogana(37)  who  divided  the  region  of 
integration  into  a large  number  of  rectangles  and 
evaluated  the  integral  by  quadrature  at  each  step 
of  an  iteration  process. 


Analogous  results  for  the  same  airfoil  in  a 
slightly  supersonic  stream  with  = 1.15  are 
shown  in  Fig.  5(10,27).  The  calculations  were 
performed  using  the  two  relaxation  methods  and  also 
a finite-difference  solution  of  the  Euler  Eqs . (5 ) by 


Fig.  5.-  Pressure  distributions  and  bow  wave  loca- 
tions for  a thin  circular-arc  airfoil  indicated 
by  FCR  and  NCR  solutions  of  the  transonic  small 
disturbance  equations,  and  an  F.uler  solution. 


Corresponding  results  for  a NACA  0012  airfoil 
are  presented  in  Fig.  7 together  wit(i  an  essentially 


Fig.  7.-  Pressure  distributions  for  an  NACA  0012 
airfoil  indicated  by  integral  equation  and  FCR 
solutions  of  the  transonic  small  disturbance 
equations,  and  by  a numerical  solution  of 
the  Euler  equations. 
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exact  solution  of  the  Euler  Eqs.  ( U ) C 52 ) for  the 
subcritical  case.  The  x-1  /2  singularity  in  u^ 
that  would  he  present  if  the  usual  thin-airfoil 
boundary  condition  Eq.  (8)  were  employed  has  been 
avoided  in  the  integral  equation  results  of  Fig.  7 
by  using  Eq.  (lU)  for  the  boundary  condition. 

Aside  from  modest  differences  near  the  shock,  the 
integral  equation  results  are  virtually  identical 
to  those  of  the  finite  difference  relaxation 
method  and  obtained  with  considerably  less 
computation. 

Attention  is  turned  now  to  unsteady  flow. 

Fig.  8 presents  pressure  distributions  for  an 
NACA  6l*Al|10  airfoil  oscillating  in  pitch  about  a 
mean  angle  of  attack  of  2°  in  a flow  with 


The  results  Just  described  are  of  great  theo- 
retical significance;  however,  their  computing  cost 
is  too  high  for  the  engineering  analysis  of  aero- 
elastic  and  flutter  problems  which  requires  a large 
number  of  cases  to  be  considered  individually  since 
solutions  cannot  be  superposed.  Simplification 
may  be  achieved  by  turning  to  the  small  disturbance 
theory  for  unsteady  transonic  flow,  but  the  need  to 
consider  every  case  separately  remains  since  Eq.(7) 
is  nonlinear.  For  many  applications,  it  is  fruit- 
ful to  decompose  the  problem  into  steady  and 
unsteady  components  by  letting 

Z = Z + Z,  C = C + C (15) 
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Fig.  8.-  Time-dependent  solution  of  Euler  Eq.  (5) 
for  the  pressure  distribution  on  an  NACA  6Aa2i10 
airfoil  oscillating  in  pitch. 


= 0.72,  as  calculated  by  Magnus  and  Yoshihara 
(21)  using  a time-dependent  method  to  solve  the 
Euler  Eqs.  (3).  For  the  three  smallest  reduced 
frequencies  k = UC/2U,,,,  the  amplitude  is  2°; 
but  it  is  only  0.4°  for  k = 2.5.  For  each  fre- 
quency, the  pressure  distributions  on  the  upper 
and  lower  surfaces  are  shown  for  the  minimum  and 
maximum  angles , and  for  the  angle  of  attack  both 
increasing  and  decreasing  through  the  mean.  At 
even  the  smallest  frequency,  k =0.1,  the  results 
differ  distinctly  from  those  of  quasi-stationary 
theory  (k  = 0).  As  k increases,  the  amplitude 
of  the  pressure  variations  diminishes,  and  the 
shock  movement  virtually  vanishes  for  k =2.5. 
The  latter  may  be  of’  importance  to  those  devel- 
oping approximate  theories,  since  it  suggests 
that  shock  movement  may  be  disregarded  for  high 
frequency  oscillations. 


where  the  barred  quantities  satisfy  Eqs.  (7) 
through  (9)  with  the  time  derivatives  omitted.  The 
remaining  unsteady  quantities  should  satisfy 
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for  the  boundary  conditions  and  pressure  relation 
given  by  Eqs.  (8)  and  (9)  with  tildes  over  <$>,  Z, 
and  Cp,  but  significant  simplification  can  be 
achieved  for  small  amplitude  motions  of  an  airfoil 
with  nonvanishing^steady-state  disturbance  field 
by  disregarding  4>x4*xx  to  obtain  a linear  differ- 
ential equation  for  Although  the  dependence 

of  the  variable  coefficients  $x  arid  on  the 

non linear  steady-state  solution  makes  the  equations 
more  difficult  to  solve  than  those  of  ordinary 
linearized  theory,  Cp  is  linearly  dependent  on  the 
amplitude  of  the  unsteady  motion,  and  results  for 
various  unsteady  motions  about  a single  steady- 
state  condition  can  be  superposed  to  determine 
solutions  for  more  complicated  motions. 

Fig.  9 presents  results  calculated  in  this  way 
by  fc’hlers(53)  using  a finite  difference  relaxation 
method  for  an  NACA  614A006  airfoil  at  zero  angle  of 
attack  with  an  oscillating  flap  extending  over  the 
rear  quarter  chord.  The  results  are  qualitatively 
similar  to  the  experimental  data(5M,  but  the  un- 
steady pressures  are  generally  overpredicted. 
Traci(55)  et  al. , have  confirmed  Ehler's  results 
by  application  of  essentially  the  same  method  to 
the  low  frequency  counterpart  of  Eq.  (lb)  in  which 
is  omitted  in  addition  to  the  nonlinear  term. 
Perhaps  the  most  notable  feature  of  the  results  is 
the  unsteady  pressure  peak  at  about  midchord  in  the 
supercritical  example  with  H»  = 0.853  that  is 
completely  unheralded  in  the  resuits  for  the  sub- 
critical  example  with  Mac  = '.79**.  The  peak  at 
three-quarter  chord  is  that  normally  expected  at 
the  hinge  line  for  any  subsonic  flow.  The  second 
peak  in  the  supercritical  case  is  associated  with 
the  accumulation  of  receding  waves  as  they  move 
upstream  from  the  flap  until  they  are  slowed  and 
extensively  refracted  by  the  zone  of  supersonic 
flow.  The  data  for  = 0.90  and  0.96  indicate 
that  such  effects  disappear  as  expected  when  the 
flow  upstream  of  thi-  hinge  line  becomes  supersonic. 

Supplementing  the  development  described  above 
is  a search  for  more  economical  approximate  net  nodi. 
Notable  among  these  is  an  extension'  ^,5  ( • C{*  the 
local  linearization  method  developed  original  ly  ( 39) 


20 


Local  Linearization 


Steady  state 


Steady  state 


Quasi- steady 


Fig.  10.-  Load  distribution  on  a thin  circular- 
arc  airfoil  oscillating  in  pitch  about  the  nose 


Steady  state 


axisymmetric  flow.  Although  some  of  the  analyses 
have  been  based  on  the  more  complete  theories, 
extensions  to  nonaxisymmetric  slender  bodies  and 
wings  have  been  confined  almost  exclusively  to  the 
small  disturbance  theory  described  by 


together  with  appropriate  boundary  conditions 
a pressure  relation  modified  to  include  <fy  + 
within  the  parentheses  of  Eq.  (9). 


Figure  11  shows  pressure  distributions  for  a 
10-percent  thick  parabolic-arc  body  of  revolution 
for  several  from  0.9  to  1.1  calculated  using 

the  local  linearization^^)  and  numerical  NCr(62) 
methods  together  with  experimental  data'°3).  These 
results , and  also  the  previously  unpublished  FCR 
results  calculated  at  NEAR,  are  in  good  agreement 


Pressure  distributions  on  an  NACA  6UAO06 
airfoil  with  oscillating  flap. 


for  steady  flow.  Results  for  the  magnitude, 

|dCp|,  and  phase,  $2  of  the  unsteady  aerodynamic 
loading  at  = 1 on  a 6-percent  thick  circular- 
arc  airfoil  oscillating  in  pitch  about  the  nose 
with  ampltidue  5 and  various  reduced  frequencies 
k are  shown  in  Fig.  10.  Also  included  are  the 
results  of  quasi-stationary  theory  and  linearized 
theory  to  which  the  local  linearization  results 
converge  for  small  and  large  k_.  Rather  similar 
results  have  been  obtained  by  Isogaif58)j  Kimble 
and  Wu(59)5  and  Dowelli^O)  through  application  of 
somewhat  different  approximations  based  on  the 
general  notions  of  the  local  linearization  method. 


Stint} 


■ ting  L 

- NCR 

-Local  linearization 
rcR 


11.-  Pressure  distributions  for  a parabolic 
arc  body  of  revolution  indicated  by  various 
solutions  of  the  small  disturbance 
equations,  and  by  experiment. 


Several  of  the  methods  developed  originally 
for  two-dimensional  flow  have  been  extended  to 
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over  most  of  the  body,  but  notable  differences 
exist  near  the  rear.  Some  of  these  may  be  attri- 
buted to  viscous  phenomena,  but  effects  of  the 
cylindrical  sting  model  support  and  wind  tunnel 
walls  are  also  significant.  To  demonstrate  the 
effects  of  the  sting,  NCR  results (62)  are  presented 
for  a shape  that  conforms  to  the  test  model  and 
sting  combination.  The  results  are  in  improved 
agreement  with  the  data,  but  significant  differ- 
ences remain. 

To  investigate  the  effects  of  the  walls,  Bailey 
(62)  carried  out  NCR  calculations  for  the  body- 
sting combination  in  a circular  wind  tunnel  having 
the  same  cross-section  area  relative  to  the  body 
as  the  Ames  1^-foot  wind  tunnel  where  the  tests 
were  conducted.  To  simulate  the  ventilated  walls 
of  the  test  section,  calculations  were  performed 
for  a porous  wall  with  various  porosities.  The 
results  for  a porosity  of  0.5  are  shown  in  Fig.  12 
together  with  those  for  free  air  and  for  an  open 


Fig.  12.-  Pressure  distributions  and  surface 
pressure  drag  for  a parabolic-arc  body  of 
revolution  with  sting  in  free  air,  an  open 
jet,  and  a wind  tunnel  with  porous  walls. 

jet.  Sedin  and  Karlsson^1^)  have  obtained  similar 
results  for  the  porous  wall  using  their  alternat- 
ing direction  integration  method.  The  theory 
indicates  that  the  walls  produce  an  upstream  shift 
of  the  shock  wave,  thereby  bringing  the  calculated 
pressure  distribution  and  pressure  drag  into  vir- 
tually perfect  agreement  with  the  measurements. 
Prudence  should  be  exercised  in  accepting  these 
results  as  definitive,  however,  since  the  boundary 
condition  applied  at  the  wall  is  a highly  simpli- 
fied representation  of  a complex  local  flow. 

It  has  long  been  known (3)  that  the  solution 
for  certain  thickness-dominated  steady  transonic 
flows  past  slender  wings,  bodies,  and  wing-body 
combinations  can  be  decomposed  into  four  simpler 
components  by  use  of  the  transonic  equivalence 
rule.  As  illustrated  by  Fig.  13  with  4>xt  and 
f(y*z**x,t)  set  to  zero,  three  of  these  require 
solution  of  only  the  two-dimensi Dnal  Laplace  equa- 
tion for  flows  in  each  transverse  y,z  plane;  and 
the  fourth  requires  solution  of  the  transonic 
small  disturbance  equation  for  axisymmetric  flow 
past  an  ’’equivalent”  body  of  revolution  having  the 
same  longitudinal  distribution  of  cross-section 
area  as  the  original  aerodynamic  shape.  Cheng  and 
Hafez  ^5)  and  Barnwell (66)  have  extended  that 
result  to  lift-dominated  cases,  for  which  it  is 
necessary  to  include  additional  higher-order 
(multipole)  inner  solutions  beyond  the  first-order 
thickness  (source)  and  lift  (dipole)  solutions. 


Fig.  13.-  Schematic  representation  of  unsteady 
transonic  equivalence  rule  for  slender  wing-body 

combinations  for  mi idly  unsteady  motions. 

The  first-order  thickness  and  lift  inner  solutions 
satisfy  the  two-dimensional  Laplace  equation;  but 
the  higher-order  solutions  satisfy  Poisson’s  equa- 
tion in  which  the  right-hand  side  f(y,z;x)  is  a 
function  of  the  lower-order  solutions.  We  have 
recently  extended  the  analysis  to  unsteady  flows (57 )f 
arriving  at  the  result  portrayed  in  f ig.  13. 

Figure  lU  presents  results  of  an  application  of 
the  equivalence  rule  for  steady  flow  with  Mac  = 1 
to  slender  bodies  with  elliptic  cross  sections  and 
the  same  longitudinal  distribution  of  cross-section 
area  as  a parabolic-arc  body  of  revolution  of  thick- 
ness ratio,  1/1 2(°7'.  The  solution  for  the  equi- 
valent body  was  determined  using  the  local 


Fig.  lU.-  Pressure  distributions  on  two  parabolic- 
arc  bodies  of  elliptic  cross  section;  = 1. 

linearization  method,  the  only  method  available  at. 
the  time  these  results  were  originally  published, 
and  are  for  free-air  flow  past  a complete  body, 
i.e.,  without  a sting.  After  making  allowances 
for  deficiencies  over  the  rear  of  the  body  that 
can  be  attributed  to  effects  of  wind  tunnel  walls 
and  model  support,  the  results  show  that  the  equi- 
valence rule  provides  a simple  and  accurate  means 
for  treating  transonic  flow  past  slender  bodies  of 
noncircular  cross  section. 

Figure  15  presents  a similar  comparison  of 
theoretical  and  experimental  results  for  Ma>  = 1 
for  one  of  the  elliptic  cross-section  bodies  of 
Fig.  lU  at  angles  of  attack  of  2°,  U 0 , and  6°'67)# 
Good  agreement  is  again  found  along  most  of  the 
body,  with  notable  differences  occurring  on  the 
rear.  The  latter  could  be  rectified  by  replacing 
the  solution  for  the  equivalent  body  by  a solution 
like  that  of  Fig.  IP  that  accounts  for  effects  of 
the  model  support  and  wind-tunnel  walls. 

In  contrast  to  the  foregoing,  the  development 
of  methods  for  solving  the  nor  linear  equations  for 
unsteady  axisymmetric  and  three-dimensional  tran- 
sonic flows  has  barely  begun.  We  have  applied  the 
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Fig.  15.-  Pressure  distributions  on  the  body 
of  Fig.  lU  with  a/b  =2  at  three  angles 
of  attack,  a:  M =1. 

local  linearization  method  to  flow  with  = 1 
past  a slender  body  of  revolution  engaged  in  a 
variety  of  unsteady  motions.  Fig.  16  reproduces 

results(57,68)  por  the  amplitude  and  phase  of  the 
surface  pressures  on  a nonlifting  body  undergoing 


shown  to  agree  with  local  linearization  solutions 
for  steady  flow  past  three  such  bodies.  Extension 
of  the  local  linearization  analysis  to  the  pitching 
motion  of  slender  bodies  of  revolution  at  Mo,  = 1 
has  been  carried  out  recently (69 ) and  used  to  deter- 
mine unsteady  surface  pressures  and  static  and 
dynamic  stability  derivatives  for  conical  and  para- 
bolic-arc bodies.  The  results  are  presented  to- 
gether with  those  of  other  approximate  theories  and 
are  shown  to  agree  well  with  available  data. 

IV.  Three-Dimensional  Flows  Past  Wings 
and  Wing-Body  Combinations 


Much  can  be  learned  by  the  study  of  two-dimen- 
sional and  slender-body  flows,  but  rational  aircraft 
design  requires  solutions  for  three-dimensional 
aerodynamic  configurations.  Fortunately,  the  finite 
difference  relaxation  method  can  be  generalized  to 
these  situations,  and  modern  computers  are  able  to 
perform  the  computations  in  an  acceptable  time  for 
steady  flows  past  certain  classes  of  three-dimen- 
sional shapes.  To  date,  the  analysis  has  been 
based  almost  exclusively  on  the  small  disturbance 
theory,  which  in  its  simplest  form  is  described 
by  Eq.  17  with  <t>xi  and  <f>tt  omitted. 

Figures  17  and  18  present  comparisons  of  calcu- 
lated and  experimental  data  for  a swept  ONERA  M6 
wing  at  an  angle  of  attack  of  3°  for  M =0.81) 
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Fig.  17--  Pressure  distributions  on  a swept  ONERA 
M6  wi ng ; M = 0 . 8U , a=  3°. 


Fig.  16.-  Local  linearization  solution  for 
pressures  on  the  forepart  of  a parabolic-arc 
body  of  revolution  in  oscillatory  dilatation 
proportional  to  the  local  body  radius; 

T = 0.10,  = 1. 

small  oscillatory  dilatations  proportional  to  the 
local  body  radius.  As  in  the  two-dimensional 
applications,  the  results  approach  those  of  lin- 
earized theory  for  high  frequency  oscillations, 
and  quasi-stationary  theory  for  low  frequencies. 

The  latter  is  illustrated  in  the  plot  in  the  upper 
right  in  which  the  pressure  distributions  indicated 
by  the  unsteady  analysis  for  a basic  body  of  thick- 
ness ratio,  T = 0.10,  undergoing  slow  expansion 
to  r = 0.12  and  contraction  to  t = 0.08  are 


and  0.93(70).  Computations  were  performed  using 
both  the  NCR  and  FCR  methods  and  required  about 
5 CPU  minutes  on  the  CDC  7600  computer.  As  in  the 
two-dimensional  applications,  the  FCR  method  pre- 
dicts a more  downstream  location  for  the  shock  wave 
than  the  NCR  method,  with  the  difference  being  more 
marked  at  the  higher  Mach  number.  Again,  the  data 
agree  better  with  the  NCR  calculations,  but  this 
is  certainly  due  to  a fortuitous  cancellation  of 
errors  in  the  shock  jumps  by  disregarded  viscous 
effects. 

The  experimental  data  also  reveal  a weak  oblique 
shock  upstream  of  the  main  shock  and  swept  back 
about  35°,  but  no  indication  of  it  is  given  by  the 
calculations.  This  is  to  be  expected  because 
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As  an  alternative  to  the  finite  difference 
approach  to  three-dimensional  transonic  flow  anal- 
ysis, mention  should  be  made  of  the  alternating 
direction  integration  method  of  Sedin  and  Karlsson 
(6^),  which  has  recently  been  shown  to  be  capable 
of  providing  plausible  pressure  distributions  on 
wings  and  wing-body  combinations.  Further  investi 
gation  of  the  accuracy  and  computational  require- 
ments of  the  method  is  necessary  to  determine  its 
competitiveness  with  the  finite  difference  methods 


V . Helicopter  Rotors 


The  onset  of  transonic  flow  over  the  outer  por- 
tion of  the  rotor  blades  is  one  of  the  primary 
conditions  that  sets  the  performance  limits  of 
modern  helicopters.  A transonic  regime  near  a 
blade  tip  is  inherently  nonlinear,  three-dimensional 
and,  with  forward  flight,  unsteady.  As  illustrated 
in  Fig.  2C)(9),  the  shock  on  the  upper  surface  moves 
downstream  at  points  B and  C,  is  nearly  stationary 
at  D,  and  moves  upstream  at  E.  Subsequently , the 
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18.-  Pressure  distributions  on  a swept 
ONERA  M 6 wing;  = 0.93,  a = 3°. 


Eq.  (17)  provides  a poor  approximation  for  shocks 
with  sweep  angles  greater  than  about  20°.  A sub- 
stantial improvement  can  be  made,  however,  bv 
adding  two  terms  to  Eq.  (17)  to  obtain 71/ 


Quasi-steady 

Unsteady 


for  steady  transonic  flow.  Studies  of  this  and 
related  extensions  of  the  small  disturbance  theory 
are  currently  in  progress.  Moving  closer  to  a 
complete  airplane  configuration.  Fig.  19  presents 
experimental  and  NCR  calculated  pressures  for  a 


Fig.  20.-  Pressure  distributions  for  a simulated 
two-dimensional  helicopter  rotor. 


shock  propagates  off  the  front  of  the  airfoil  into 
the  oncoming  flow.  In  the  calculation  of  the  theo- 
retical results  of  Fig.  20,  the  unsteady  motion  of 
the  rotor  blade  in  forward  flight  was  simulated  by 
an  NACA  0012  airfoil  undergoing  simultaneous  sinu- 
soidal Mach  number  and  incidence  variations  with  a 
phase  difference  of  l80°.  The  governing  equation 
is  a two-dimensional  low  frequency  version  of 
Eq.  (17)  with  ijyj  and  omitted  and  with 

Mq.  + M^sin  i(/,  which  represents  the  instantaneous 
free-stream  Mach  number,  in  place  of  M^  in  the 
coefficients  of  (t^  and  iji^  and  the  boundary 
conditions.  The  quantity  Mp  is  the  tip  Mach 
number  due  to  rotation , and = Sit  is  the  indicated 
azimuth  angle  on  the  rotor  disk.  The  quasi-steady 
results  were  computed  by  solving  the  same  equation 
with  <(>xt  = 0,  so  that  time  appeared  only  as  a 
parameter.  The  shaded  region  indicates  the  portion 
of  the  cycle  for  which  the  flow  is  subcritical,  and 
for  which  the  unsteady  values  for  the  pressures  are 
nearly  quasi-steady.  In  the  supercritical  region, 
however,  the  pressures  and  shock  wave  location  lag 
the  quasi-steady  results . 


Fig.  19.-  Pressure  distribution  on  a wing-body 
combination;  M = 0.93,  a = 0°. 


swept-wing- fuselage  configuration  at  M^  = 0.93 
and  a = 0°(70).  The  agreement  with  experiment  is 
good  along  the  fuselage  centerline  and  the  two  in- 
board stations.  In  the  computed  results,  the  wing 
root  shock  propagates  laterally  to  y/s  = 0.60, 
but  the  experimental  shock  dissipates  before 
reaching  this  point.  The  source  of  the  disagree- 
ment is  not  clear  but  may  be  due  to  viscous 
effects  or  to  deficiencies  associated  with  swept 
shock  waves. 
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The  results  of  Fig.  20  correspond  to  strip 
theory,  but  three-dimensional  effects  are  obviously 
important  for  a helicorter  rotor,  particularly  near 
the  tip.  Caradonna  and  Isom(72)  have  investigated 
this  feature  of  rotor  aerodynamics,  and  calculated 
unsteady  FCK  solutions  reproduced  in  Fig.  21  for 
both  two-  and  three-dimensional  representations. 


Fig.  21.-  Two-  and  three-dimensional  calculations 
of  pressure  distributions  for  several  ij;  at  97-7 
percent  of  the  rotor  radius  with  blades  of  aspect 
ratio  = 15  and  6-percent  thick  circular-arc  pro- 
files operating  with  an  advance  ratio,  U^/fJR  = 0.U, 
and  a tip  Mach  number,  = 0.93,  at  i|i  = 90°. 

Comparison  shows  that  the  three-dimensional  solu- 
tion indicates  that  the  flow  is  far  less  expanded 
(the  peak  negative  Cp  is  less  in  the  supersonic 
region  and  the  shock  wave  does  not  reach  the 
trailing  edge),  the  return  to  subsonic  flow  occurs 
earlier,  and  no  upstream  propagating  wave  is  indi- 
cated for  this  case.  Caradonna  and  Isom  state 
that  such  a wave  is  seen  in  their  results  for  a 
higher  Mach  number  three-dimensional  calculation, 
but  its  amplitude  is  about  an  order  of  magnitude 
less  than  indicated  by  the  two-dimensional  calcu- 
lation. This  brief  account  serves  to  indicate 
that  helicopter  aerodynamics  provides  a new  field 
of  application  for  transonic  flow  theory.  The 
available  results  are  very  incomplete,  but  they 
illustrate  the  significance  of  transonic  flow  for 
helicopter  rotors,  and  the  need  for  a more 
comprehensive  study. 

VI.  Rotating  Turbomachinery 

Flow  through  a high-speed  fan  or  compressor  is 
three  dimensional,  can  include  complex  shock 
systems,  and  is  unsteady  even  in  a rotating  frame 
of  reference  if  a complete  stage  consisting  of  a 
rotor  and  stator  or  a fan  preceded  by  inlet  guide 
vanes  is  considered.  Effects  of  viscosity  and 
turbulence  are  important , particularly  in  the  aft 
stages,  and  impact  of  the  wake  of  one  stage  on  the 
blades  of  the  next  is  an  important  source  of  vibra- 
tion and  noise.  The  calculation  of  transonic  flow 
through  high-speed  turbomachinery  must  be  one  of 
the  most  formidable  problems  in  aerodynamics,  but 
the  technological  need  exists  and  methods  are 
beginning  to  be  developed. 

In  current  work,  the  traditional  approach  is 
taken  in  which  the  inviscid  flow  field  is  consid- 
ered first,  leaving  the  viscous  and  turbulence 
effects  to  be  added  later,  hopefully  as  small  or 
localized  perturbations.  Since  variations  of 
pressure  and  flow  direction  may  be  substantial, 
the  small  disturbance  theory  is  not  expected  to 


be  as  useful  as  it  is  for  the  external  aerodynamics 
of  the  airplane.  The  appropriate  equations  for 
the  inviscid  analysis  are,  therefore,  the  Euler 
equations,  or  possibly  the  complete  potential 
equation  if  the  shock  waves  are  not  too  strong. 

In  the  most  advanced  analysis,  Erdos  et  al.,'73) 
have  approximated  the  three-dimensional  problem 
by  a pair  of  two-dimensional  problems,  which  are 
intended  ultimately  to  be  interacted  to  obtain  the 
complete  solution.  Figure  22  presents  the  essen- 
tial ideas  of  this  decomposition  and  of  the 


Fig.  22.-  Coordinate  systems  and  grid  networks 
for  interacting  two-dimensional  flows  used  to 
simulate  three-dimensional  flow  through  a 
rotor-stator  stage  of  a turbomachine. 

coordinate  system  evolved  from  it.  The  hub-to- 
casing  solution  accounts  for  effects  of  blockage 
of  the  flow  area  by  the  blades  and  the  anticipated 
boundary  layer  as  well  as  the  geometry  of  the  hub 
and  casing,  but  not  for  any  circumferential  varia- 
tions of  the  flow.  It  provides  the  flow  properties 
in  a meridional  plane  as  averaged  from  one  blade  to 
the  next,  the  coordinates  R = r(m)  of  the  curved 
axisymmetric  stream  surface  upon  which  the  blade- 
to-blade  analysis  is  carried  out,  and  thereby  the 
spacing  b(m)  between  two  arbitrarily  selected 
adjacent  stream  surfaces.  The  problem  that  must 
be  solved  for  the  blade-to-blade  surface  represents, 
therefore,  a two-dimensional  flow  past  a series  of 
blades  in  the  curvilinear  coordinate  system  (m,9) 
indicated  in  Fig.  22  with  the  radius,  R(m),  and 
thickness,  b(m),  of  the  stream  sheet  provided  by 
the  solution  of  the  hub-to-casing  solution. 

Expressed  in  a coordinate  system  rotating  with 
the  blades  at  an  angular  velocity  V.  so  as  to 
obtain  a steady  flow  problem  if  interaction  with 
other  blade  rows  is  disregarded,  the  equations  to 
be  solved  are 


3p  3(pV)  3(pU)  _ PU  d(Rb) 

3t  + R30  + 3m  - Rb 

3(pV)  3(p  + pV2 ) 3(pUV)  ..  PUV  d(Rb) 

3t  R36  3m  Rb  dm 

-f  (V  + 2ftR)  £ 

3(PV)  . 3(PW)  , 3(ptpUz)  _ _ PUZ  d(Rb) 

3t  R30  3m  Rb  dm 


+ p(VtliR)2  dR 
R dm 
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3(pEr)  3(pVHp)  3(pUHr)  PUHr  d(Rb) 

3t  + R30  + 3m  Rb  dm 


(19) 


where  U and  V represent  velocity  components  in 
the  m and  0 directions,  and  in  which  the  rela- 
tive total  energy  and  total  enthalpy  are  defined 
by 


E = E - W2R,  H = H - VfiR  (20) 

r r 

The  right-hand  sides  of  Eq.  (19),  which  differ 
from  those  of  Erdos  et  al.(73)  because  of  correc- 
tion of  errors  in  the  transformation  to  curvi- 
linear  coordinates,  arise  from  streamwise  varia- 
tions of  thickness  and  radius  of  the  stream  sheet; 
they  vanish  if  the  equations  are  applied  to  a two- 
dimensional  cascade  flow.  The  corresponding 
equations  for  a nonrotating  inertial  system  may 
be  obtained  by  setting  ft  to  zero. 

Only  a few  results  obtained  by  solving  these 
equations  have  been  reported,  and  they  are  rather 
provisional  because  of  the  relatively  small  number 
of  grid  points  used  in  the  finite  difference 
calculations.  Results  of  blade-to-blade  calcu- 
lations of  Erdos  et  al.,(73)  are  presented  in 
Fig.  23  together  with  experimental  data  for  a 
high-speed  (1500-fps)  fan  tip  section  obtained 


Fig.  23.-  Experimental  and  calculated  pressure 
distributions  for  flow  through  a 1500-fps  rotor. 

from  an  array  of  fast  response  pressure  gages 
mounted  on  the  casing  wall.  The  test  fan  was 
preceded  by  a set  of  guide  vanes  and  followed  by 
a row  of  stators,  b\;t  the  unsteady  interactions 
were  neglected  in  the  calculations  and  only  the 
effects  of  the  rotor  was  considered.  Although 
the  grid  network  was  very  coarse,  the  calculated 
flow  bears  a recognizable  resemblance  to  the  ob- 
servations, particularly  the  indications  of  an 
oblique  shock  from  the  leading  edge  of  the  upper 
blade  that  reflects  off  the  lower  blade  and  reim- 
pinges on  the  upper  blade  near  the  trailing  edge. 
These  developments  give  promise  of  substantial 
improvements  in  predictive  capabilities  for  the 
design  and  analysis  of  high-speed  turbomachinery. 
The  flows  are  very  complex,  however,  and  continued 
effort  will  be  required  for  some  time  before 
accurate  and  economical  methods  will  be  available 
for  the  routine  solution  of  these  problems. 


VII.  Concluding  Remarks 

The  preceding  discussion  has  provided  a review 
of  recent  developments  in  steady  and  unsteady 
transonic  aerodynamics.  Many  references  are  cited; 
but  many  significant  contributions  are  not,  or  are 
hidden  in  references  to  summary  papers.  In  many 
instances,  the  selection  has  been  made  on  the  basis 
of  interlocking  relations  that  help  to  evaluate  the 
various  results,  and  to  provide  a continuing  base 
from  which  further  discussion  can  proceed  smoothly. 

Overall,  it  is  clear  that  tremendous  advances 
are  being  made  in  the  analysis  of  transonic  flows, 
and  that  these  problems  no  longer  appear  as  for- 
midable as  they  once  did.  Indeed,  some  of  the 
simpler  two-dimensional  and  axi symmetric  steady 
flows  may  be  considered  solved,  with  even  alter- 
native methods  available  to  choose  among.  The 
research  frontier  is  moving  now  to  more  complex 
steady  and  unsteady  three-dimensional  flows  past 
wings,  wing-body  combinations,  helicopter  rotors, 
and  through  turbomachinery  fans  and  compressors. 

The  modem  computer  has  brought  immense  calculating 
power  to  bear  on  these  problems,  and  the  goal  of 
replacing  the  wind  tunnel  with  a computer  is 
beginning  to  look  more  achievable  than  ever  before. 
However,  all  will  not  fall  into  place  by  itself. 
Much  work  must  be  done,  but  the  directions  are 
indicated  and  the  rewards  of  improved  aerodynamic 
design  and  analysis  are  sufficient  to  demand  the 
effort  be  made. 
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